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Intense  Relativistic  Electron  Beam  Interaction 
with  a Cool  Theta  Pinch  Plasma 


I.  INTRODUCTION 

The  recent  availability  of  electron  beam  generators1  5 capable  of 
producing  in  excess  of  1 MJ  of  relativistic  electrons  in  times  of  ~ 100 
ns  has  stimulated  a great  deal  of  interest  in  their  application  to  Con- 
trolled Thermonuclear  Fusion  Research.  The  beam-plasma  systems  under 
investigation  can  be  divided  into  inertially  confined  and  magnetically 
confined  configurations.  In  the  inertial  confinement  case,  present  ef- 
forts are  directed  toward  the  development  of  beam  generators  capable  of 
producing  the  required  high  power  (~  1014  w),  short  pulse  (~  10  s) 
electron  beams,  and  the  focusing  of  such  beams  onto  small  targets.  7 
On  the  other  hand,  in  magnetically  confined  systems,  it  is  total  energy, 
rather  than  power,  that  matters,  so  long  as  the  beam  pulse  duration  is 
less  than  the  energy  loss  time  of  the  confinement  system.  Therefore, 
the  10  s pulsed  intense  relativistic  electron  beams  (iREB's^  pres- 
ently under  development^  are  appropriate,  and  it  appears  that  the  main 
problem  areas  here  lie  in  the  compatibility  of  IREB  heating  and  the  con- 
finement system.  In  particular,  the  strength  of  the  interaction  must  be 
such  as  to  efficiently  deposit  the  beam  energy  in  the  plasma.  In  an 
open  system,  such  as  a high  magnetic  field,  long  solenoid,"  the  energy 
deposition  length  must  be  the  system  length.  If  a beam  can  be  injected 
Note:  Manuscript  submitted  December  27,  1976. 
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into  a toroidal  system10  without  seriously  affecting  the  energy  confine- 
ment, the  deposition  length  can  be  much  longer  so  long  as  the  beam 
energy  is  deposited  before  it  is  lost  by  such  processes  as  synchrotron 

radiation.11  Therefore,  the  strength,  as  measured  by  the  energy  deposi- 

. 

tion  length  (or  time!,  and  the  characteristics  of  the  interaction  between 
an  intense  beam  and  a plasma  will  determine  the  potential  applicability 
of  IREB 1 s to  heating  magnetically  confined  plasmas.  In  particular,  since 
the  classical  Maxwellian  plasma  interaction  lengths  are  too  long  for  effi- 
cient deposition  in  plasmas  even  in  the  1017  cm  density  range,  collec- 
tive energy  coupling  processes  are  required. 

Beam-plasma  interaction  experiments  previously  reported1  20  have 
observed  interaction  strengths  which  imply  that  a collective  interaction 
must  be  taking  place.  The  present  work  was  directed  toward  investigating 
the  beam-plasma  interaction  under  conditions  in  which  the  strength  and 
the  characteristics  of  the  physical  processes  involved  could  be  stidied 
in  detail.  It  was  designed  so  that  several  of  the  difficulties  in 

I 

interpretating  the  results  of  previous  experiments  were  eliminated.  Thus, 
a 4 m long  uniform  plasma  was  used  to  avoid  magnetic  field  and  plasma 
inhomogeneities.  Secondly,  the  plasma  density  was  high  enough  that 
Thomson  scattering  could  be  used  to  determine  plasma  electron  density  and 
temperature.  Finally,  the  beam  current  density  was  kept  low  enough 
2 kA/cnf  ) that  detailed  local  magnetic  field  measurements  could  be 
made  within  the  beam  channel  without  probe  destruction. 

Summarizing  the  experiment,  a 0.5  - 1 MeV,  25  - '0  kA,  oO  - 70  ns 
electron  beam  was  injected  into  preionized  plasma  confined  in  a 4 m long 
theta  pinch  by  a 4-5  kC.  magnetic  field.  In  addition  to  Thomson  scattering 
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and  magnetic  probes,  diagnostics  used  to  study  the  beam  plasma  interac- 
tion included  1)  diamagnetic  loops  surrounding  the  plasma,  principally 
to  compare  with  the  other  diagnostics,  2)  visible  light,  to  diagnose  the 
characteristics  of  the  plasma  discharge,  3)  x-band  microwave  apparatus, 
to  monitor  radiation  near  the  electron  cyclotron  frequency,  and  4)  hard 
x-rays,  to  study  the  angular  spread  of  the  beam.  Two  different  plasma 
conditions  were  investigated.  In  the  first  one,  100  mTorr  of  helium  was 
partially  ionized  to  produce  a (0.5  - 4)  x 1 o14/  cm3,  1 - 2 eV  plasma  at 
the  time  of  beam  injection.  For  these  experiments,  the  beam  to  plasma 
density  ratio  was  ~ 10  - 10  and  the  fraction  of  beam  energy  coupled 

to  the  plasma  was  ~ 5'fVm.  However,  a large  fraction  of  this  energy  was 
lost  to  atomic  processes  (ionization  and  line  radiation).  Preliminary 
results  for  these  plasma  conditions  were  presented  by  Goldenbaum,et  al ,24 
The  second  plasma  condition  was  highly  ionized  hydrogen  at  a density  of 
(0.5  - 4)  x 10x“/cm'  and  temperature  of  2 - 3 eV  (electrons  and  ions). 

In  this  case  the  beam  to  plasma  density  ratio  was  10  - 10  1 and  the 

energy  coupling  efficiency  was  1 - 2$/m.  Since  electron  and  ion  temper- 
atures were  equal  at  the  time  of  beam  injection,  the  ion  acoustic  in- 
stability was  probably  not  present.  The  high  plasma  density  and  small 
neutral  fraction  made  it  possible  to  obtain  the  heating  rate  during  the 
beam  pulse  by  Thomson  scattering.  For  plasma  density  below  2 x 10'  “/cm' , 
evidence  for  nonclassical  heating  was  obtained.  Furthermore,  detailed 
magnetic  probe  measurements  were  made  of  changes  in  magnetic  field  com- 
ponents during  the  beam-plasma  interaction.  The  axial  field  showed  a 
rapidly  rising  (~  20  ns)  diamagnetic  signal  which  (on  average)  contin- 
ued to  rise  slowly  throughout  the  beam  pulse.  The  heated  cross  sectional 
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area  after  beam  passage  was  more  than  twice  the  beam  area.  The  azimuthal 
magnetic  field  indicated  a net  axial  current  density  within  the  inter- 
action region  which  was  much  higher  than  predicted  by  return  current 
theory  assuming  classical  dissipation.27  The  anomalously  large  heating 
rate  and  the  high  net  current  density  can  be  explained  by  approximately 
the  same  effective  collision  frequency  during  the  beam  plasma  interac- 
tion. Finally  at  the  highest  densities  (>  2 x lO^/cm3),  classical  re- 
sistivity return  current  heating  adequately  accounts  for  the  observations. 
Part  of  the  results  obtained  under  the  highly  ionized  plasma  conditions 
were  reported  in  preliminary  form  by  Dove,  et  al.25 

Clearly  an  understanding  of  these  results  requires  a comparison 
with  expectations  based  upon  the  various  beam  plasma  interaction  mecha- 
nisms which  have  been  discussed.26  We  find  that  the  electron-electron 
two  stream  instability  can  provide  the  anomalous  resistivities  and  heat- 
ing rates  observed  in  the  present  experiments.  If  the  model  used  to  ex- 
plain our  results  applies  into  the  high  temperature  plasma  regime,  then 
it  can  be  expected  that  longer  pulse  duration  electron  beams  will  have 
greater  overall  plasma  heating  efficiency.  This  is  because,  for  a given 
beam- to-plasma  density  ratio,  the  instantaneous  coupling  efficiency  in- 
creases as  the  plasma  is  heated. 

The  organization  of  this  article  is  as  follows:  In  order  to 

facilitate  comparison  with  our  results,  in  Sec.  II  we  review  intense 
beam-plasma  interaction  mechanisms,  and  previous  experimental  results. 

In  Sec.  Ill,  we  describe  the  experimental  apparatus,  including  the  plas- 
ma source,  the  electron  beam  generator,  the  diagnostics  and  the  charac- 
teristics of  the  initial  plasma  and  the  beam.  In  Sec.  IV,  we  discuss 
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the  experimental  results,  and  in  Sec.  V they  are  interpreted  in  terms  of 
theoretical  predictions.  Finally  in  Sec.  VI,  we  discuss  the  implications 
of  the  present  work  to  application  of  IREB's  to  controlled  fusion  in 
magnetically  confined  systems. 

II.  REVIEW  OF  INTERACTION  MECHANISMS  AND  PREVIOUS  EXPERIMENTS 
A . Theoretical  Interaction  Mechanisms 

Coupling  of  the  energy  of  an  IREB  into  a plasma  by  collective 
processes  can  occur  by  several  mechanisms.  Several  of  these  have  been 
reviewed  recently  by  Breizman  and  Ryutov.‘?°  We  divide  the  collective 
mechanisms  into  two  main  categories,  microscopic  and  macroscopic.  In  the 
first,  the  beam  excites  an  instability,  and  the  individual  electrons 
interact  directly  with  large  amplitude  waves  which  are  present  at  satura- 
tion. Thus,  the  beam  transfers  its  energy  to  waves,  which  pass  it  in 
turn  to  the  plasma.  The  instability  which  has  received  the  most  atten- 
tion in  this  regard  is  the  electron-electron  two  stream  instability,2''-35 
although  other  instabilities  have  also  been  considered 6»3S  The  macro- 
scopic category  includes  effects  which  are  unique  to  intense  beams, 
namely  the  induced  return  current, 27,3 7-40  the  presence  of  large  self 
fields  if  the  beam  is  not  fully  charged  neutralized  or  current  neutral- 
ized by  plasma  motion, 26*37t41 »42  and  the  large  transverse  pressure  ex- 
erted by  the  beam  against  a confining  magnetic  field,  particularly  if  the 
beam  is  rotating  across  field  lines .2j ,43 »44  We  now  proceed  to  discuss 
briefly  some  of  the  main  characteristics  of  these  mechanisms. 

The  electron-electron  two  stream  instability  (e-e  mode)  is 
driven  by  the  relative  drift  between  beam  and  plasma  electrons.  Linear 


5 


growth  rates  for  unstable  waves  have  been  calculated,"  recalculated,”1 
reviewed,^'  and  then  calculated  again'  ‘ ’3o  for  evermore  "realistic"  con- 
ditions, meaning  conditions  which  more  closely  approximate  the  experi- 
ments. For  example,  the  maximum  growth  rate  6 for  a "cold"  beam  satis- 
fying ( n^/n  ) 3 y « 1 is 


3 cu  / it 
8 - \ - — - 


2 V \ 2 n 


(cos1'  ? + sin2  §)  , 


where  n^  and  n^  are  the  beam  and  plasma  densities,  respectively,  vine" 
is  the  beam  electron  energy,  ou  = (n^e^/e^m)  is  the  plasma  frequency 
and  $ is  the  angle  of  the  wavevector  k with  respect  to  the  beam  propaga- 
tion direction.  (The  quantities  m,  -e,  c,  and  eQ  are  the  mass  and  charge 
of  the  electron,  the  velocity  of  light  and  the  permittivity  of  free 
space,  respectively.)  A beam  is  "cold"  when24 


Av  /c  =?  § + AE/^mc2  « 


where  the  parallel  velocity  spread  Av  is  due  to  either  the  beam  energy 

II 

spread  AE  or  the  spread  in  angle  of  beam  electron  velocity  vectors  rela- 
tive to  the  beam  propagation  direction,  represented  here  by  the  mean 
angle  9.  This  is  the  so-called  hydrodynamic  limit  of  the  instability, 
in  which  all  of  the  beam  particles  see  the  same  phase  of  the  unstable 
waves,  and  interact  with  them  coherently  as  one  fluid,  on  the  timescale 
of  the  instability  (~  6 X).  In  the  other  limit,  called  the  kinetic,  or 
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warm  beam,  limit,  the  growth  rate  is  given  by 


.28 


% 1 


£ _ rj  b_  p 

P vn  9"  at2  + k2c& 


(3) 


where  k = | k | sin  § and  it  is  assumed  that  0 » AE/v^mc.2.  Similar 

growth  rates  can  be  obtained  for  strong  beams  ^(n^/n^)1^3  y » ']  in 
these  two  limits,'*  ‘ and  for  the  case  when  a magnetic  field  is  present, 
in  which  case  other  instabilities  are  also  expected.'  Unfortunately, 
the  linear  growth  rate  tells  only  a part  of  the  story  since  it  is  the 
nonlinear  limit  to  which  the  instability  goes  which  determine  the  inter- 
action strength  and  characteristics.  Several  limiting  mechanisms  have 
been  discussed,  such  as  quasilinear'  > " > ' and  several  wave-wave  scat- 
tering processes'1  ’ and  each  gives  its  own  interaction  length  > 

The  energy  can  be  removed  from  the  beam  either  directly  via  beam  elec- 
tron-wave interaction,  as  demonstrated  graphically  by  computer  simula- 
tions, 0 j34-3  or  by  generating  anomalous  resistivity,33  or  both.  Al- 
though the  predicted  interaction  lengths  vary  by  many  orders  of  magnitude 
they  are  all  short  compared  to  classical  interaction  length  unless  the 
plasma  is  too  inhomogeneous  in  the  direction  of  beam  propagat ion.""  For 
example,  the  quasilinear  length,'*31  which  is  the  shortest  one,  is 


c n 

10  — -E 


kT 


at  n^  me 


- ve2 

2 


(M 


where  kT  is  the  plasma  electron  temperature.  For  a lO^/cm"  density, 
e 

few  eV  plasma  and  a 101'  /cm’1  beam  with  any  mean  angle  0,  It  is  small 
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compared  to  c/uj^,  an  extremely  short  length.  In  fact  it  is  unreasonably 
short  since  the  wave  energy  density  used  to  derive  it  is  sufficiently 
large  that  the  weak  turbulence  approximations  used  in  the  derivation  are 
invalid  ,3~’  We  note  that  for  most  of  the  experiments  performed  to 
date,  the  kinetic  growth  rate,  Eq.  (3),  should  apply,  although  Thode34 
has  obtained  reasonably  good  correlation  of  theory  with  experimental  re- 
sults with  a modified  version  of  the  hydrodynamic  formulation. 

While  the  microscopic  processes  just  discussed  depend  upon  the 
ratio  n^/n  and  not  on  total  beam  current  1^,  the  macroscopic  interac- 
tions may  have  1^  dependence.  Therefore,  it  is  useful  to  define  the 
ratio  v/y>  by  which  the  strength  of  an  IREB  is  commonly  measured,  before 
proceeding  further.  The  magnitude  of  the  current  in  an  electron  beam  is 
given  by 


4tt  e me3  Nec 


I,  h NeV,  = 
b b 


^TTe  me 
o 


i3b  ■ 


(5) 


where  N is  the  number  of  electrons  per  meter  of  beam  length  and  p,  = V,  /c. 

b b 

In  terms  of  the  quantity  v,  the  number  of  electrons  per  classical  elec- 
tron radius  of  beam  length  e2/4rreomc2  =2.83  X 10-1_  m,  this  is 


I,  = 17000  vP,  amperes, 
b b 


05) 


As  a measure  of  the  strength  of  the  beam  self  fields,  we  determine  the 

current,  I , at  which  the  cyclotron  radius  r of  an  electron  at  the 
c c 

boundary  of  a cylindrical  beam  of  radius  b equals  b/2.  At  this  current, 

the  electron  is  strongly  affected  by  the  beam's  self  magnetic  field  Ba(b). 

0 
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Since  r = ymV, /eB  (b)  and  B (b)  = [a  I /2rrb  = I /2nbe  c2,  we  find 
c b 9 S ob  b o 

I = 2ne  c2bB„(b) 
c o 9 ' 

4rre  me3 

= 2 PbY  = 17000  Pby  • (7) 

e 

This  "critical"  current41’40  means  that  when  v/y^  1>  beam  self  fields, 

if  unneutralized,  severely  affect  electron  motion,  whereas  for  v/y  « 1 

they  do  not.  Typical  value  of  I are  28  kA  at  500  keV  kinetic  energy 

(y=2),  50kAat  1 MeV  (y  = 3)  and  85  kA  at  2 MeV  (y  = 5). 

Turning  now  to  the  macroscopic  interaction  processes,  we 

first  discuss  the  induced  return  current.  It  is  a result  of  the  fact 

that  intense  beams  are  pulsed.  The  dl/dt  in  the  beam  front  produces  a 

dB/dt  which,  by  Maxwell's  equation  V x E = - 5B/dt,  gives  rise  to  an 

electric  field.  (The  electric  field  and  magnetic  induction  are  E and  B, 

respectively.)  This  field  tends  to  slow  down  beam  electrons,  but  it 

also  tends  to  accelerate  plasma  electrons  back  down  the  beam  channel  so 

as  to  eliminate  the  3B/3t  (Lenz  Law).  For  a collisionless  plasma,  when 

« n^  and  utpb/Vb  » 1,  current,  neutralization  is  virtually  complete 

throughout  the  beam  cross  section.  If  oj  b/V.  < 1,  then  neutralization 

P b 

of  the  beam  current  density  within  the  beam  channel  is  only  partial. 

These  results  are  correct  in  the  absence  of  applied  magnetic  f ields,"'7' " >40 
or  in  the  presence  of  longitudinal  or  transverse  fields,"  so  long  as 

uJl  > of,  where  ui^  = eBQ/m  is  the  cyclotron  frequency  of  plasma  electrons 
in  the  applied  field  Bo< 
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T , , where  t , is  the  momentum 
col  col 


Since  a plasma  is  not  collisionless,  the  induced  plasma 

current  diffuses  out  of  the  beam  channel.  This  has  been  shown"7  to 

occur  on  a timescale  T = b'  /Cc/uu  ) | T ,,  where  t , is  the  momentum 

p J col’  col 

transfer  collision  frequency  of  the  plasma  electrons.  T » in  a 

fully  current  neutralized  beam  (cu  b/v^  » 1)  occurs  because  energy  is 
delivered  inductively  from  beam  electrons  to  plasma  electrons  to  make  up 
for  the  energy  dissipated  (converted  to  plasma  thermal  energy)  by  Ohmic 

heating. 26,31 t4G>47  t0  summarize,  the  energy  deposited  per  unit  length 

47 

in  the  plasma  by  the  return  current,  Q,  is  given  by 


t b 

■H : 


2TTr^  • Edr  = / dtRI^  , 


where  I is  the  plasma  current,  is  the  beam  current  density,  and  R is 
the  plasma  column  resistance  per  unit  length  within  the  beam  channel. 
Cylindrical  symmetry  and  axial  uniformity  are  assumed.  If  there  are  no 
space  charge  electric  fields  present,  then  Q is  the  difference  between 
the  work  done  by  the  beam  and  the  magnetic  field  energy  per  unit  length 
at  t ime  t 


- f L H dt  - i 
J 


LI?(t)  , 


where  I=I^+l(=0att=0  assuming  initially  complete  current  neu- 
tralization), and  L is  the  inductance  per  unit  length  of  the  entire  beam 
plasma  system.  Since  1^  « V^,  and  Y « ( 1 - V^/c2)  1 implies 
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6v./V,  « by/y , where  6v,  is  the  velocity  decrease  of  a beam  electron 
b b b 

due  to  an  energy  change  of  by,  a relativistic  beam  electron  can  loose  a 
large  fraction  of  its  energy  without  changing  the  beam  current  signifi- 
cantly. Therefore,  we  may  take  1^  out  of  the  integral  in  Eq . (9)  to 
obtain 


Q = LI,  I - - LI f I.  - - I 

b 2 lb  2 


Then  the  energy  deposited  by  the  end  of  the  beam  pulse  is  determined  by 
the  net  current  in  the  system  at  that  time.  A reasonable  limiting  value 
for  the  net  current  is  I since  a larger  current  is  going  to  face  severe 
magnetic  self  pinching  as  we  have  seen.  The  time  required  for  the  cur- 
rent in  a beam  with  I,  > I to  reach  I = I is  ~ T( I /I,  ),  at  which  time, 

be  c c b ’ 

Q « Ll|(v/y  - 1/2) . This  is  always  a shorter  time  than  for  beam  energy 
deposition  by  binary  collisions  between  beam  electrons  and  plasma  parti- 
cles for  beams  and  plasmas  of  interest.  How  much  shorter  depends  upon 
Tco^,  whether  it  is  due  to  coulomb  collisions  or  is  enhanced  by  the 
presence  of  microturbulence  due  to  instabilities.  This  microturbulence 
must  be  low  frequency  if  it  is  to  affect  plasma  electrons,  and  can  occur 
either  parametrically  as  a result  of  the  e-e  mode''  ’3  or  due  to  insta- 
bilities generated  by  the  relative  drift  between  plasma  electrons  and 
ions  resulting  from  the  return  current  flow.  In  the  latter  category  are 
such  instabilities  as  electron- ion  two  stream  and  ion  acoustic.'  > 1>4  ,47 


The  interaction  length  from  this  process  can  be  estimated  as 
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(y  - l)mc2  (y  - l)mc2T 


(1 


T dI 

eL  — 


eLl( t ) 

P 


dt 


where  is  the  beam  pulse  duration.  It  will  vary  directly  with  the  ef- 
fective collision  time  tco^  since  (to  first  approximation)  I will  vary 
inversely  with  tcq^.  Note  that  the  energy  stored  in  the  magnetic  field 
1/2  LI2  in  Eq.  ( 10)  also  ultimately  ends  up  in  the  plasma  by  Ohmic  dissi 
pation.26’ * However,  this  process  can  take  much  longer  time  than  the 
dissipation  during  the  beam  pulse  since  instabilities  have  presumably 
become  much  weaker  or  even  stabilized. 

Also  included  in  the  macroscopic  collective  interaction  pro- 
cesses category  were  the  effects  of  the  beam's  large  self  fields  and 
transverse  pressure.  Suppose,  for  example,  the  beam  is  injected  into  a 
low  pressure  neutral  gas.  Then  until  the  gas  is  turned  into  a plasma 
the  beam  self  electric  and  magnetic  fields  can  build  up.  This  can 
strongly  affect  the  characteristics  of  the  gas  breakdown  process,  deter- 
mining, in  turn,  the  beam  and  plasma  characteristics  after  breakdown  anc 
the  type  of  interaction  that  can  take  place.  These  processes  can 
reasonably  be  expected  to  be  strongly  affected  by  the  beam  strength 
(v/y)*  Tw°  examples  have  been  studied.  In  one,  the  beam  self  magnetic 
field,  in  the  absence  of  an  applied  guide  field,  can  cause  a strong 
pinch,  ■’  heating  ions  as  well  as  electrons.  On  the  other  hand,  if  ai 
IREB  is  injected  into  a plasma  or  neutral  gas  with  a substantial  frac- 
tion of  its  energy  in  the  transverse  direction,  its  transverse  pressure 
can  substantially  exceed  the  confining  pressure  of  the  applied  field. 
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This  nonequilibrium  situation  can  result  in  the  generation  of  large 
amplitude  magnetosonic  waves  which  can  also  heat  both  ions  and  elec- 
trons.43 Such  expansion  waves  can  also  be  driven  by  hot  plasma  if,  for 
example,  the  plasma  electrons  are  strongly  heated  by  some  instability 
mechanism  (e.g.  - electron-electron  two  stream)  in  a time  short  compared 
to  the  characteristic  plasma  expansion  time.''1  All  of  these  mechanisms 
have  been  observed  in  experiments,  as  we  shall  see  in  the  experimental 
review. 

B . Previous  Experimental  Results 

The  availability  of  intense  electron  beams  in  the  mid  1960's 
soon  resulted  in  experimental  studies  in  which  beam  propagation  character- 
istics were  investigated  in  plasmas  and  in  initially  neutral  gasesl1 * 56 
However,  no  attempt  was  made  to  study  the  beam- to- plasma  energy  coupling. 
The  first  experiment  designed  to  do  this  was  that  of  Altyntsev,  et  al.lg 
Although  a relatively  weak  beam  was  used,  substantial  nonclassical  beam- 
to-plasma  energy  coupling  was  observed  for  n^/n^  in  the  range  1 to  10 
In  fact  in  all  studies  of  energy  deposition  by  an  IREB  in  a plasma  in 
which  the  density  ratio  was  ~ 10  or  higher,  including  in  the  experiment 
reported  here,  energy  deposition  greater  than  can  be  explained  by  classi- 
cal processes  has  been  observed.  Moreover,  even  when  the  density  ratio 
was  small  (again  including  the  present  experiment),  classical  dissipation 
of  the  beam- induced  return  current  appears  to  be  the  energy  coupling 
mechanism.  We  note  that  even  this  is  a collective  ( nonturbulent , macro- 
scopic) interaction  mechanism  that  would  not  be  present  in  an  ordinary 
beam-plasma  interaction. 
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Returning  now  to  the  experiment  of  Altyntsev,  et  al . , 
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2-5  MeV,  10  kA,  50  ns  electron  beam  was  injected  into  a 5 m long, 
magnetic  mirror  confined  hydrogen  or  argon  afterglow  plasma  in  the  densi- 
ty range  1011  - 1014/cm3.  The  beam  density  was  approximately  1011/cm 
in  the  interaction  volume,  and  the  midplane  magnetic  field  was  'S  2.5  kG. 
The  interaction  strength  was  diagnosed  by  monitoring  the  beam  propagation 
efficiency  to  the  end  of  the  interaction  region  by  beam  calorimetry,  and 
by  using  diamagnetic  loops  to  measure  the  total  transverse  energy  per 
unit  length  of  the  particles  (plasma  electrons  and  ions,  and  beam  elec- 
trons). The  interaction  was  found  to  be  strongest  when  the  beam  and 
plasma  density  were  comparable,  with  most  of  the  beam  energy  not  reaching 
the  end  of  the  system.  However,  the  diamagnetic  loops  indicated  that 
10$  of  the  d irected  beam  energy  was  converted  into  transverse  particle 
energy,  this  quantity  being  optimized  at  a plasma  density  of  ~ 1C1,  cm 
Where  the  rest  of  the  energy  went  was  unknown.  Thus,  the  interaction 
length  for  beam  loss  was  -3m  for  a large  enough  beam-plasma  density 
ratio.  For  coupling  of  beam  energy  into  the  plasma,  it  was  perhaps  10  m, 
much  longer  than  the  quasilinear  length  given  in  Eq . (4) , but  still 
orders  of  magnitude  shorter  than  is  possible  by  classical  collisional 
processes.  Note  that  10$  of  the  beam  energy  equally  distributed  among 
all  plasma  particles  in  the  system  corresponds  to  tens  of  keV  per  elec- 
tron ion  pair. 

In  the  Altyntsev,  et  al . , experiment  the  beam  was  relatively 
weak,  with  v/y  < 0.1.  Similar  experiments  were  performed  by  Smith1  and 
Okamura,  et  al.,"  ' using  even  lower  v/y  beams,  but  again  observing  much 
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stronger  than  classical  beam-plasma  interaction  (also  using  beam  calorim- 
etry and  diamagnetic  loops).  Abrashitov,  et  al.,1  and  Arzhannikov,  ^t 
al.,17  followed  up  the  work  of  Altyntsev,  et  al . , with  a more  completely 
diagnosed  experiment,  but  still  a v/v  « 1 beam.  In  particular,  at  a 
plasma  density  ~ 1014/cm3,  Thomson  scattering  was  used  to  obtain  the 
plasma  transverse  electron  temperature  kT^.  The  resulting  n^kT^  was  ~ 5 
times  smaller  than  the  plasma  transverse  energy  inferred  from  diamagnetic 
loops.  Other  diagnostics  in  this  and  other  experiments  suggest  a hot 
electron  component,1  J’ 2 or  hot  ions,1  or  both,  may  account  for  this 
difference.  Since  diamagnetic  loop  measurements  were  made  within  100  ns 
of  beam  injection,  other  possible  explanations  for  the  discrepancy  are 
the  residual  effects  of  beam  diamagnetism,  and  magnetosonic  waves  which 
have  not  yet  damped.44’31  The  largest  signals  were  observed  for  n^/n^  = 
0.05  - 0.005.  However,  at  the  lower  value  of  n ( n,  /n  > 0.05),  beam 

Pop 

propagation  was  poor  (Fig.  4b,  Ref.  17). 

Two  experiments  followed  shortly  after  the  work  of  Altyntsev, 
et  al . , in  which  stronger  beams,  v/ v « 1 - 5»  were  used.  Kapetanakos  and 
Hammer13  injected  a 20  - 40  kA,  400  keV,  50  ns  beam  into  a 40  cm  long, 
mirror  confined  afterglow  hydrogen  or  helium  plasma  in  the  density  range 
1011  - 1014/cm3.  Miller  and  Kuswa2  J used  a 50  kA,  350  keV  beam  of  30  ns 
duration,  and  a 30  cm  long  1012  - 1014/cm  plasma  confined  in  a uniform 
magnetic  field.  Both  experiments  had  beam  densities  of  101'  /cm  , both 
used  diamagnetic  loops  as  the  principal  diagnostic  for  beam  energy  trans- 
ferred to  the  plasma,  and  both  observed  maximum  values  of  transfer  ef- 
ficiency of  ~ 5$.  However,  these  maxima  were  at  different  density 
ratios  in  the  two  experiments:  at  n^/np  1 for  Kapetanakos  and  Hammer, 
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and  at  n^/n^  ~ 10  ‘ for  Miller  and  Kuswa.  This  difference  may  have  been 
due  to  the  different  beam  characteristics  or  plasma  confinement  configura- 
tions. Kapetanakos  and  Hammer  also  found  plasma  diamagnetism  to  be  in- 
dependent of  magnetic  field  above  2.5  kG  and  to  scale  as  B'  below  that 
field,  probably  a beam  or  plasma  confinement  effect.  Miller  and  Kuswa 
observed  soft  x-rays  at  the  lower  plasma  densities  confirming  the  presence 
of  substantially  heated  plasma  electrons.  They  also  made  the  first  men- 
tion of  magnetosonic  oscillations.’1 

Further  experiments  at  intermediate  values  of  v/ v were  per- 
formed by  Korn,  e t a 1 . » and  Ekdahl,  et  al.j  on  the  same  apparatus.  A 
10  - SO  kA,  550  keV,  j0  ns  beam  (v/va:  1/2  - 3)  was  injected  into  a 
fully  ionized  plasma  in  the  density  range  lO1^  - 5 X 101/cm  . The  mid- 
plane magnetic  field  in  this  1.  m long  experiment  was  typically  2.7  kG. 

As  in  the  case  of  Kapetanakos  and  Hammer,  the  coupling  efficiency  was 
highest,  up  to  ~ 20f*  of  the  beam  enetjy  deposited  in  the  plasma,  at  the 
higher  beam  to  plasma  density  ratios,  again  using  diamagnetic  loops. 
However,  a neutral  particle  detector  was  also  used  and  it  was  found  that 
the  plasma  ions  had  gained  a substantial  amount  of  energy.  This  was 
probably  due  to  large  amplitude  magnetosonic  waves  since  diamagnetic  loop 
oscillations  scaled  as  B/n  as  predicted  by  theory.'  ’ In  one  set  of 
experiments,  the  energy  deposition  as  a function  of  current  was  con- 
sistent with  return  current  heating  in  the  presence  of  ion  sound  turbu- 
lent resistivity.1 " In  the  other  experiments,2  the  heating  rate  re- 
quired 10  times  that  resistivity.  The  use  of  a foilless  diode13  “ in  the 
first  set,  of  experiments  and  an  ordinary  foil  diode  in  the  second  may 
be  the  explanation  for  the  different  interaction  characteristics  here. 
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suggested  that  most  of  the  energy  deposition 


Furthermore  Korn,  et  al . ,1 c 
occurred  in  the  first  few  tens  of  cm  of  plasma.  "Previous  workers  had 
also  seen  evidence  that  this  might  be  the  case. 12> 14, la, it  yet  ^pub- 
lished work  on  the  same  apparatus  by  Sethian,  et  al.,°9  has  shown  a very 
strong  dependence  of  the  interaction  strength  upon  beam  mean  angle  as 
determined  by  scattering  in  the  anode  foil.59  This  result,  obtained  by 
both  Thomson  scattering  and  diamagnetic  loops  at  a plasma  density  of 
5 X 1013/cm3,  points  to  the  presence  of  the  electron- electron  two  stream 
instability,34  although  return  current  heating  may  still  be  present  as 
well.50 

Turning  now  to  a higher  v/y  beam  experiment,  Miller20  injected 
a 600  kA,  100  ns  beam  (v/y  r*  5)  into  a ~ 2 x 1013/cm3  density  plasma. 

Two  parameters  were  varied  - the  anode  foil  thickness  and  the  neutral  gas 
pressure.  The  first  of  these  determines  the  beam  electron  mean  angle 
relative  to  the  propagation  direction,  an  important  parameter  if  the  two 
stream  instability  is  operative . 2 “>33 »34  A decrease  by  a factor  of  ~ 10 
was  observed  in  diamagnetic  loop  signals  for  the  thicker  foils.  The  in- 
jected beam  current  density  was  decreased  from  a maximum  of  ~ 6 kA/cm^ 

(%  « 1012/  cm3)  by  a factor  of  ^ 2 by  the  thicker  foils.  Another  im- 
portant result  here  was  that  beam  diamagnetism,  certainly  substantially 
increased  by  the  thicker  foils,  was  not  dominating  the  diamagnetic  loop 
signals.  The  change  in  neutral  gas  density  holding  the  plasma  density 
constant  should  not  affect  plasma  heating  mechanisms  unless  collisions 
with  neutrals  become  so  numerous  as  to  damp  instabilities,  or  the  density 
changes  significantly  during  the  pulse.  The  small  decrease  of  the  ob- 
served diamagnetic  loop  signal  at  the  higher  pressures  may  have  been  due 
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pulse,  but  also  ionization  and  line  radiation  can  result  in  substantial 
plasma  kinetic  energy  loss  until  full  ionization  is  achieved.  The 
maximum  diamagnetic  loop  signal  implied  a 1 - 25=  energy  deposition  effi- 
ciency, about  l/k  to  1/2  that  expected  based  upon  earlier  experi- 
ments1^’14 19  at  the  same  value  of  n^/n^.  Note  that  we  expect  the  high 
v/y  beam  Miller  used  to  have  no  particular  effect  on  interaction  strength 
since  n^  ^ 20  n^  and  the  potentially  large  beam  self  fields  will  be 
easily  neutralized  hy  plasma  motion. 

In  initially  neutral  gas  experiments,  however,  the  high  v/y 
can  be  important.  In  the  case  of  VanDevender,  et  al . ,^1  a 550  keV,  kA 
(v/y  = 1.6)  100  ns  beam  was  focused  to  25  kA/cm^  (n^  « 5 X 1012)  and  in- 
jected into  hydrogen  in  the  pressure  range  50  mTorr  - 10  Torr  with  no 
applied  magnetic  fields.  At  the  lower  pressures,  rather  substantial  net 

BQ , corresponding  to  a net  current  of  ~ 506  of  the  primarv  beam  current, 

a 

was  observed.  Since  the  net  current  is  in  the  direction  of  the  beam 
current,  during  the  beam  pulse  the  plasma  current  is  opposite  to  the  net 
current.  Therefore,  the  plasma  is  antipinched  by  the  ^ X IB  volume  force. 
After  the  beam,  the  plasma  current  reverses  and  the  j y B force  is  in- 
ward. This  antipinching  followed  by  pinching  was  observed  by  streak 
photography,  and  an  ion  temperature  of  ~ 100  eV  was  found  by  visible 
spectroscopy.  The  electrons,  however,  were  found  (by  Thomson  scattering) 
to  have  a non-Maxwe Ilian  distribution  function  with  mean  energy  only 
~ 10  eV  and  a density  ~ lO^/cm3.  When  the  plasma  was  preionized,  the 
interaction  was  much  weaker,  and  results  were  consistent  with  return 
current  heating,  both  with  and  without  anomalous  resistivity  depending 
upon  plasma  conditions. 
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Prono,  et  al . ,gg  did  use  an  applied  magnetic  guide  field  in 
their  high  v/y  neutral  gas  experiment.  They  injected  a 1 MeV,  200  - 500 
kA  beam  (v/y  = 4 - 10)  of  150  ns  duration  into  neutral  hydrogen  ranging 
from  JO  mTorr  to  1 Torr.  In  the  range  30  - 130  mTorr,  where  gas  break- 
down occurs  slowly,42  the  13  kA/cm2  beam  (n^  « 2.5  X 1012/cm3)  inter- 
acted very  strongly,  losing  ~ 50^  of  its  energy  in  the  1 m long  experi- 
ment, and  producing  diamagnetic  loops  signals  corresponding  to  ~ 1 j/cirr  . 
Even  a small  amount  of  preionization  substantially  reduced  the  interac- 
tion, pointing  to  the  importance  of  the  initially  neutral  gas.  Similar- 
ly, operation  at  the  higher  pressures,  where  gas  breakdown  occurs  quick- 
ly,'* also  greatly  weakened  the  interaction. 

The  final  previously  reported  beam-plasma  interaction  experi- 
ment to  be  discussed  here,  by  Kapetanakos,  et  a l.,*3  utilized  a sharp 
magnetic  cusp  to  nonadiabatically  convert  a beam  propagating  parallel  to 
a magnetic  field  (500  keV,  20  - 40  kA,  50  ns)  to  a beam  rotating  across 
magnetic  field  lines.  This  rotating  beam  of  density  ~ 5 X 1011/cm  was 
then  injected  into  a plasma  in  the  density  range  10 12  - 1014/cm‘.  Con- 
trary to  the  parallel  propagation  case,  this  experiment  found  the  strong- 
est beam- to-plasma  energy  coupling,  again  as  determined  by  diamagnetic 
loops,  to  be  at  a plasma  density  of  ~ 1014/cm \ This,  however,  was  in 
agreement  with  predictions  by  Chu  and  Rostoker.  1 In  their  model,  the 
plasma  heating  is  caused  by  the  fields  associated  with  the  cross  field 
return  current.  Ion  heating  was  indicated,  although  not  measured,  in 
this  experiment  due  to  the  strong  magnetosonic  oscillations  which  were 
observed  above  101  /cm  plasma  density. ** 
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We  have  seen  that  all  of  the  beam-plasma  interaction  experi- 
ments reported  so  far  have  observed  rates  of  energy  deposition  in  the 
plasma  which  imply  the  presence  of  one  or  more  collective  mechanisms. 

The  results  of  some  of  the  experiments1'2-1  7,13  ~al  indicated  the  presence 
of  the  electron-electron  two  stream  instability,  particularly  at  the 
lower  plasma  densities,  and  others  seemed  to  imply  return  current 

heating.  In  all  of  the  experiments  both  mechanisms  could  have  been  pre- 
sent. One  difficulty  in  interpreting  most  of  the  experimental  results 
is  the  lack  of  detailed  local  measurements  of  plasma  conditions  after 
beam  injection.  Of  the  two  experiments  reported  in  which  Thomson  scat- 
tering was  used  to  obtain  plasma  electron  density  and  temperature,  one 
was  performed  in  initially  neutral  gas,23  and  the  other  was  in  a plasma 
fn  which  conditions  would  allow  both  electron-electron  two  stream  and  ion 
acoustic  instabilities  to  be  present.17  The  remaining  experiments  de- 
pended upon  nonlocal  measurements,  mainly  diamagnetic  loops,  to  infer 
beam  plasma  coupling,  and  if  this  is  done  shortly  after  beam  injection, 
beam  diamagnetism  and  magnetosonic  waves  can  contribute  significantly  to 
the  signals. 

III.  EXPERIMENTAL  APPARATUS  AND  CHARACTERIZATION  OF  THE  BEAM  AND  INITIAL 
PLASMA 

The  electron  beam  generator  used  in  the  present  experiment  is  an 
8 Q,  water  dielectric  coaxial  Blumlein5  pulse  forming  line  driven  by  a 
Marx  generator.3  The  latter  consists  of  twelve  0.5  uF  capacitors  which 
were  charged  as  high  as  80  kV  each,  at  which  voltage  a total  of  19  kJ  is 
stored.  When  the  Marx  generator  is  switched  into  a series  configuration, 
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an  output  gap  Is  closed  and  the  pulse  forming  line  is  charged  to  peak 
voltage  (up  to  1.2  MeV)  in  approximately  1 |js . At  this  time  four 
pressured  gas  switches  are  simultaneously  triggered  by  an  externally  pro- 
duced high  voltage  trigger  pulse,  and  a 60  ns  ( fwhm)  pulse  is  produced. 

The  command  triggering  of  the  pulse  line  switches  was  incorporated  to 
allow  precise  timing  between  beam  injection  into  the  plasma  and  the 
Thomson  scattering  diagnostic  (see  below).  The  multiple  switching  serves 
to  reduce  the  effective  switch  inductance  and,  therefore,  the  rise  time 
of  the  output  pulse.  The  negative  voltage  pulse  from  the  pulse  forming 
line  is  applied  to  a standard  vacuum  field  emission  diode,  of  the  type 
described  by  Parker,  et  al.6g  In  this  experiment  cathodes  consisted  of 
flat  carbon  discs  4 - 7 */2  cm  in  diameter.  The  anode  was  typically  a 
titanium  foil  25  um  thick,  and  was  spaced  0.7  - 1.2  cm  from  the  cathode. 

Operating  voltages  on  the  cathode-anode  gap  ranged  from  500  keV  to  as 
high  as  1.4  MeV.  (Note  that  a Blumlein  pulse  forming  line  produces  a 
pulsed  voltage  equal  to  the  charging  voltage  into  a matched  load,  8 Q in 

this  case.  Twice  the  charging  voltage  is  produced  on  an  open  circuit.5)  i 

Electron  beam  currents  ranged  from  25  kA  to  125  kA . Further  details  of 

the  generator  are  available  elsewhere.'0  Figure  1 shows  a sample  set  of 

oscilloscope  traces  of  a.)  the  Marx  generator  charging  the  Blumlein 

pulse  forming  line,  b.)  the  electron  beam  diode  voltage,  and  c.)  the 

electron  beam  current  in  the  diode. 

Magnetic  field  coils  were  mounted  on  the  accelerator  diode  so 
that  electrons  leaving  the  cathode  followed  magnetic  field  lines  through 
the  anode  and  into  the  theta  pinch  target  plasma.  The  rise  time  of  this 
guide  field  was  long  enough  (5  ms)  to  penetrate  the  stainless  steel, 
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brass,  and  carbon  parts  which  were  inserted  within  it  (see  Fig.  2).  The 
amplitude  of  the  guide  field  was  chosen  to  match  the  field  at  the  cathode 
anode  gap  to  the  theta  pinch  magnetic  field  at  the  time  of  beam  injec- 


The  plasma  source  in  this  experiment  was  a 4 m long,  5 Us  rise 
time  theta  pinch,  together  with  a Z discharge  preionizer.  The  theta 
pinch  itself  was  a 20  cm  inside  diameter  single  turn  coil  driven  by  a 
60  kJ,  20  kv  capacitor  bank.  The  plasma  was  contained  inside  a 15  cm 

_ c 

inside  diameter  glass  vacuum  vessel  which  had  a base  pressure  of  ~ 10 
Torr.  Gas  fill  pressures  ranged  from  5 to  200  mTorr  of  hydrogen,  deu- 
terium or  helium  depending  upon  the  desired  plasma  conditions.  The 
parameters  of  the  plasma  contained  in  the  theta  pinch  at  the  moment  of 
electron  beam  injection  were  determined  by  the  preionization  sequence. 

In  order  to  produce  a partially  ionized  plasma  with  n^  =“  (0.5  - 5)  X 
1014/cm3,  the  discharge  tube  was  filled  to  a gas  pressure  of  50  - 200 
mTorr,  and  the  Z discharge,  powered  by  a 1 uf  capacitor  bank  charged  to 
30  kV,  was  fired  only  a few  us  before  the  theta  pinch.  To  produce  a 
highly  ionized  plasma  the  tube  was  filled  to  5 - 10  mTorr  and  the  Z dis- 
charge was  fired  ~ 20  us  before  the  theta  pinch.  Access  to  the  beam 
generator  required  a ~ 1 m long  drift  section  between  the  anode  foil  and 
the  entrance  to  the  theta  pinch  coil.  The  magnetic  guide  field  at  the 
cathode  anode  gap  extended  over  this  distance.  To  insure  that  this 
drift  region  was  highly  ionized  at  the  time  of  beam  injection,  the  pri- 
mary Z discharge  electrodes  were  the  anode  foil  and  electrodes  placed 
just  before  the  theta  pinch  entrance,  as  illustrated  in  Fig.  2.  Shortly 
after  breakdown  in  this  drift  region,  the  ionizing  discharge  within  the 
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theta  pinch  tube  was  struck  from  the  theta  pinch  entrance  electrodes  to 
the  Faraday  cup  - calorimeter  (see  below)  which  terminated  the  beam- 
plasma  interaction  region  about  1/2  m inside  the  end  of  the  theta  pinch 
tube.  In  some  of  the  highly  ionized  plasma  experiments,  a 3 kJ  low  in- 
ductance capacitor  bank  was  switched  into  the  theta  pinch  coil,  the 
ringing  discharge  of  which  preheated  the  plasma  for  subsequent  further 
heating  and  compression  by  the  theta  pinch  discharge.  Figure  3 presents 
a sample  oscilloscope  trace  of  the  light  from  the  plasma  in  the  dis- 
charge tube  for  the  highly  ionized  case. 

For  the  bulk  of  the  work  to  be  presented  here,  a typical  electron 
beam  pulse  launched  from  the  diode  was  2-3  kJ.  Of  this  energy  ~ 75^ 
actually  entered  the  theta  pinch,  the  loss  occurring  mainly  at  the 
transition  between  the  slow  guide  field  and  the  fast  theta  pinch  field 
regions.  Within  our  shot-to-shot  reproducibility,  ± 15^,  all  of  the 
beam  energy  injected  into  the  theta  pinch  was  collected  by  the  Faraday 
cup  - calorimeter  3 X/2  m downstream.  The  latter  diagnostic,  of  the  type 
described  by  Pellinen,'  4 allows  beam  current  and  total  beam  energy  to  be 
determined.  Surveys  with  witness  plates  at  various  axial  locations  with- 
in the  discharge  tube  showed  the  beam  to  maintain  its  initial  cross- 
sectional  area  (e.g.,  40  cm2  with  the  7^/2  cm  diameter  cathode).  For 
example  at  the  axial  position  of  the  laser  diagnostic  it  was  typically 
distorted  into  an  elliptical  form  (~  3 cm  x ~ 6 cm)  with  its  major  axis 
in  a horizontal  plane  (the  plane  of  the  theta  pinch  slot),  and  displaced 
upward  and  toward  the  theta  pinch  slot  1-2  cm. 

Studies  of  the  beam  produced  x-ray  spectrum  and  angular  depen- 
dence from  a titanium  strip  in  the  guide  field  region  indicated  that  the 
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Scattering  in  the  0.0025  cm 
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beam  had  an  angular  spread  of  ~ 1 radian.^4 
titanium  anode  can  account  for  only  a small  portion  of  this  (~  1/3  radian 
at  the  600  keV  level  used  in  those  studies).  Therefore,  we  postulate 
that  the  angular  spread  was  caused  by  beam  interaction  with  the  plasma 
near  the  anode  foil,  or  by  magnetic  field  nonuniformities  due  to,  for 
example,  joints  in  the  guide  field  coils. 

The  electron  beam  was  injected  into  the  plasma  1 us  after  peak 
current  was  achieved  in  the  theta  pinch  coil  in  order  to  allow  the  plas- 
ma to  expand  to  a reasonably  uniform  radial  profile.  This  tended  to 
maximize  the  electron  beam  energy  which  could  be  injected.  In  the  fully 
ionized  plasma  case,  this  also  allowed  the  electrons  and  ions  to  equili- 
brate. At  the  time  of  beam  injection,  electron  and  ion  temperatures 
were  determined  to  be  2 - 3 eV,  the  former  by  Thomson  scattering  and  the 
latter  by  measuring  the  Doppler  broadening  of  the  radiation.  In  the 
partially  ionized  case  the  density  was  determined  by  an  absolute  cali- 
bration of  the  continuum  intensity  using  the  model  described  by  Griem,35 
and  by  Thomson  scattering  when  the  density  was  high  enough  (>  1014/cm“). 

The  key  measurements  of  plasma  electron  density  and  temperature 
after  beam  injection  into  the  plasma  were  made  by  Thomson  scattering'3 
of  ruby  laser  light  (69^3  A).  The  laser  had  a peak  power  of  400  MW  and 
a pulse  width  of  30  ns  ( fwhm) . The  scattering  volume  was  a 1 - lx/2  cm 
long  by  3 mm  diameter  cylinder  of  plasma  centered  on  the  axis  of  the 
theta  pinch  tube  1 1/2  m from  the  entrance.  The  scattering  vector  was 
transverse  to  the  tube  axis  and  the  applied  magnetic  field  direction  so 
that  f(v  ) was  measured.  Side  arms  in  the  glass  tube  which  passed 
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through  the  theta  pinch  coil  were  used  for  the  laser  input  and  dump. 

The  scattered  light  was  collected  via  a third  side  arm  at  90°  to  the 
other  two.  Stray  light  was  reduced  by  having  a thin  sheet  of  anodized 
stainless  steel  attached  to  the  discharge  tube  wall  opposite  the  scat- 
tered light  collection  port.  The  scattered  light  was  analyzed  with  a 5 
channel  polychromator  consisting  of  a 1/2  m spectrometer,  23  A spectral 
width  fibre  bundles  and  (RCA  7265)  photomultiplier  tubes.  The  entire 
analyzing  system  was  heavily  lead  shielded  to  eliminate  signals  due  to 
hard  x-rays  from  the  electron  beam.  Optical  and  electrical  design 
parameters  of  the  scattering  system  are  available  elsewhere.'^  The 
relative  sensitivities  of  the  scattering  channels  were  determined  using 
a Tungsten  lamp  having  a known  profile,  and  with  the  continuum  from  the 
plasma  before  and  after  beam  injection.  The  absolute  calibration  was 
then  determined  for  the  channel  on  the  laser  line  center  (6943  A)  by 
Rayleigh  scattering  from  nitrogen  gas.  Figure  4a  shows  a typical  scattered 
light  signal  on  one  of  the  channels  when  the  plasma  was  partially  ionized 
helium.  The  large  initial  signal  occurred  during  beam  injection.  It  was 
optical  continuum  during  the  beam-plasma  interaction,  not  x-rays.  This 
large  background  signal  made  density  and  temperature  measurements  by 
laser  scattering  impossible  until  about  JO  ns  after  the  start  of  beam 
injection.  Figure  4b  shows  a typical  scattering  signal  in  the  highly 
ionized  case.  The  continuum  radiation  increased  here  as  well,  but  only 
a small  amount.  Thus,  in  this  case  the  increase  served  as  a marker  for 
the  time  of  arrival  of  the  beam  at  the  laser  scattering  port.  However, 
it  did  not  prevent  density  and  temperature  measurements  by  laser  scat- 
tering during  the  beam  pulse.  Figure  5 shows  the  range  of  densities  and 
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Lemperatures  obtained  in  one  sequence  of  discharges  of  the  Z discharge- 


theta  pinch  system  in  which  the  initial  fill  was  10  mTorr  hydrogen 
(6.7  X 1014  atoms/cm3).  Although  the  density  varied  from  ~6  x 1014  cm’ 
to  ~ 4 x 1015  cm  the  temperature  remained  in  the  narrow  range  of 
2.2  - 3-2  eV.  The  large  density  variation  was  probably  due  to  nonrepro- 
ducible  desorption  of  gases  from  the  glass  tube  walls  during  the  early 
stages  of  the  theta  pinch  discharges,  and/or  nonreproducible  compression 
of  the  resultant  plasma  by  the  theta  pinch.  The  latter  could  also  ex- 
plain the  previously  mentioned  ± 15**  variation  in  beam  energy  injected 
into  the  theta  pinch. 

Visible  light  measurements  were  made  on  the  plasma  discharge  usii 
1/4  m and  1 m monochromators  together  with  phatomultipliers . For  ex- 
ample, in  partially  ionized  helium  plasmas,  the  time  history  of  He-I 
(4922  A)  and  He-II  (4686  A.)  lines  were  observed  to  determine  the  rate  of 
energy  deposition  in  the  plasma  during  the  beam  pulse.'  4 The  1 m mono- 
chromator  was  used  to  obtain  the  Doppler  broadened  (opol  A)  line 
width  in  the  highly  ionized  hydrogen  plasma.  Stark  and  Zeeman  line 
broadening  were  small  for  our  plasma  conditions,  and  resonant  charge  ex- 
change neutrals  emitted  easily  observable  levels  of  Doppler  broadened 
radiation.  The  line  width  was  measured  by  sampling  0.5  width  segment! 
of  the  line  on  successive  plasma  discharges,  three  or  more  discharges 
each  wavelength.  The  result  of  this  measurement  is  shown  in  Fig.  6. 

Also  shown  in  the  figure  is  the  instrument  profile  as  measured  separately 
with  a discharge  tube  source.  The  line  broadening,  assuming  Gaussian 

. o 

profiles,  is  about  0.54  A,  which  yields  an  ion  temperature  of  ~ 5 eV  at 
the  time  of  injection  of  the  electron  beam.  Approximately  equal  electri 
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and  ion  temperatures  is  consistent  with  collisional  relaxation  times  for 
our  conditions.  These  measurements  were  made  ~5  m from  the  entrance  to 
the  theta  pinch  coil. 

Two  different  types  of  magnetic  diagnostics  were  used  during 
various  phases  of  this  work.  Local  magnetic  field  measurements  were  made 
in  the  highly  ionized  plasma  case  with  a 3 mm  diameter,  6 turn  pickup 
loop  mounted  on  the  end  of  a 3 mm  diameter  solid  copper  outer  conductor 
50  fi  coaxial  cable.  The  loop  was  moved  vertically  across  the  plasma  on 
successive  shots  inside  a 1.25  cm  diameter,  3 tnm  wall  quartz  tube  which 
extended  all  the  way  across  the  discharge  tube  on  all  shots.  In  this  way 
the  plasma  and  beam  perturbation  by  the  probe  housing  was  the  same  on 
all  shots.  Silicon  dielectric  fluid  within  the  quartz  tube  served  to 
suppress  electrical  discharges  within  the  tube.  The  presence  of  the 
probe  housing  within  the  discharge  tube  (~  3 m from  the  entrance  to  the 
theta  pinch  coil)  was  found  to  have  negligible  effect  on  the  plasma  den- 
sity and  temperature. 

The  second  type  of  magnetic  diagnostic  device  was  a diamagnetic 
loop,  positioned  around  the  glass  discharge  chamber,  which  measures  the 
change  in  total  enclosed  axial  flux.  Because  of  the  large  voltages  in- 
duced in  the  loop  by  the  theta  pinch  discharge,  it  was  necessary  to  con- 
nect the  loop  in  series  with  a multiturn,  small  diameter  compensating 
coil  located  between  the  discharge  tube  and  the  theta  pinch  coil. 

Because  of  the  large  inductance  of  the  resulting  electrical  circuit,  the 
rise  time  of  the  loop  circuit  was  ~ 1/2  us.  Integration  time  constants 
of  5 and  10  us  were  used.  One  such  loop  was  located  ~ 1/2  m inside  each 
end  of  the  theta  pinch  coil  (see  Fig.  2). 
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IV  . 


EXPERIMENTAL  RESULTS 


A . Highly  Ionized  Case 

Several  different  experimental  runs  were  taken  injecting  the 
electron  beam  into  the  highly  ionized  hydrogen  plasma.  In  the  first  cast 
to  be  discussed,  a 60  - 75  kA  peak  current,  JO  ns  ( fwhm  current  pulse), 

~ 2 kj  beam  reached  the  Faraday  cup-calorimeter.  Cathode-anode  gap  peak 
voltage  was  1.0  MeV,  and  the  beam  cross  section  was  ~ ^+0  ciri^ . Figure  J 
shows  plasma  density  vs  temperature  data  from  Thomson  scattering  for  thi 
run.  This  figure  includes  scattering  data  taken  from  2+0  to  350  ns  after 
the  arrival  of  the  beam  front  at  the  scattering  port.  The  data  points 
shown  in  Fig.  5> obtained  by  firing  the  full  system  with  the  exception  ol 
the  electron  beam,  were  taken  interspersed  among  the  beam  shots  in  Fig. 
Therefore,  they  are  representative  of  the  plasma  into  which  the  beam  was 
injected.  Note  that  the  density  range  covered  by  the  prebeam  plasma 
data  ("prep  shots")  and  the  beam  data  is  virtually  the  same.  It  is 
clear  in  Fig.  7 that  the  temperature  is  strongly  density  dependent,  with 
the  higher  density  data  clearly  showing  the  lowest  temperatures,  and 
vice  versa,  even  though  the  scatter  is  substantial.  (The  curve  is  simp] 
to  guide  the  eye.)  Relative  errors  in  these  data  points  are  typically 
± 10$  in  both  density  and  temperature.  In  Fig.  8,  the  temperature  as  a 
function  of  time  is  plotted.  As  a result  of  the  functional  dependence 
of  temperature  on  density,  the  data  has  been  divided  into  three  density 
ranges  based  upon  the  results  shown  in  Fig.  J:  (7  ± 3)  X 1014/cm2, 

(1.5  ± -5)  X 1 01 5 / cm3  and  (3.1  ± 1.0)  X 1015/  cm3 . The  curves  drawn  for 
the  highest  density  case  will  be  discussed  in  the  next  section.  There 
appears  to  be  no  real  trend  over  the  time  period  covered.  This  is  in 
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contrast  to  our  previously  published  results2'  from  a different  series 
of  shots,  in  which  rapid  post  beam  cooling  seemed  to  occur.  Given  the 
scatter  in  the  present  data,  the  discrepancy  is  probably  due  to  the  small 
number  of  shots  in  the  previously  published  run.  The  lack  of  cooling  is 
consistent  with  the  fact  that  even  in  the  17  eV  shot,  electrons  moving  at 
the  thermal  velocity  would  take  1 Us  to  travel  half  the  length  of  the 
discharge  tube.  However,  the  present  data  does  tend  to  verify  our  pre- 
vious contention  that  most  of  the  heating  occurs  in  the  early  part  of  the 
beam  pulse.  The  data  from  the  previously  published  run  are  included  in 
Fig.  8 in  the  appropriate  density  grouping  for  comparison. 

The  series  of  shots  just  discussed  were  performed  with  the 
five  Thomson  scattering  channels  placed  symmetrically  in  wavelength 
space  around  the  line  center  at  the  exit  plane  of  the  polychromator . 

Thus,  channels  1 and  2 were  located  on  the  red  (longer  wavelength)  side 
of  69^3  A,  3 was  centered  on  69^+3,  and  channels  k-  and  5 were  on  the  blue 
(shorter  wavelength)  side.  This  allowed  a slight  shift  between  the  line 
center  in  prep  shots  and  beam  shots  to  be  observed.  Such  a shift  should 
exist  since  the  beam  leaves  behind  plasma  currents  in  the  r,  9 plane  as 
well  as  in  the  z direction  when  it  exits  the  plasma.  ' If  the  beam 
plasma  system  were  azimuthally  symmetric,  the  currents  in  the  r,  9 plane 
would  be  purely  azimuthal  and  the  plasma  electron  distribution  would 
have  a drift  in  the  theta  direction  if  the  current  is  due  to  electrons. 
Figure  9a  is  an  idealization  of  the  actual  unsymmetric  geometry  showing 
the  scattering  volume,  the  beam  position  on  a typical  shot,  and  the  dia- 
magnetic drift  direction.  Since  the  beam  cross  section  was  not  circular 
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the  system  was  even  less  symmetric  than  shown.  Figure  9b  shows  sample 
prep  shot  and  beam  shot  scattering  data  to  illustrate  the  relative 
change  between  a slight  blue  side  enhancement  in  the  prep  shot  case  to  a 
slight  red  side  enhancement  for  a beam  shot.  The  blue  side  enhancement 
in  the  case  of  prep  shots  is  a geometric  effect  determined  by  such  things 
as  the  precise  placement  of  the  fibre  bundles  in  the  exit  plane  of  the 
polychromator . Line  center,  69V5  , is  taken  to  be  at  the  average  line 

center  obtained  from  the  prep  shots  assuming  the  line  shape  is  Gaussian. 
The  calculation  is  illustrated  in  Appendix  A.  A shift  of  the  profile 
toward  the  red  side  during  beam  shots  implies  a mean  electron  velocity 
vector  <v>  such  that  k*<v>  is  negative,  where  k is  the  scattering 
vector,  the  direction  of  which  is  shown  in  Fig.  9a. 

Figure  10  shows  the  observed  line  center  for  the  beam  shots  of 
Fig.  8 as  a function  of  plasma  density  and  temperature,  together  with  the 
same  information  for  the  prep  shots  of  Fig.  5»  Errors  in  the  observed 
intensities  produce  a typical  uncertainty  in  the  deduced  line  center  in 
each  shot  of  approximately  ± 1 A.  This  is  consistent  with  the  calculated 
standard  deviation  for  the  prep  shot  line  centers  of  0.93  A.  Figures 
10a  (line  center  vs  density)  and  10b  (line  center  vs  temperature')  show 
no  apparent  density  and  temperature  dependence  for  the  line  center  posi- 
tion in  the  prep  shots.  On  the  other  hand,  for  beam  shots,  Figs.  10a 
and  10b  show  possible  trends  to  greater  line  center  for  lower  density 
and  for  higher  temperature,  respectively. 

For  the  geometric  situation  shown  in  Fig.  9a,  we  might  expect 
the  drift  velocity,  if  it  is  a result  of  plasma  electron  diamagnetic 
current  flow,  to  be  more  evident  in  the  lower  part  of  the  scattering 


volume  than  the  upper  assuming  symmetric  plasma  motion  around  the  beam 


axis . This  is  because  the  electron  drift  should  be  a minimum  near  the 
middle  of  the  beam  heated  plasma  - i.e.,  near  the  beam  axis.  Therefore, 
a series  of  beam  and  prep  shots  were  performed  looking  only  at  the  lower 
half  or  only  the  upper  half  of  the  scattering  volume.  The  results  for 
these  cases  are  shown  in  Fig.  11.  The  line  center  as  a function 
of  temperature  and  density  for  the  upper  half  of  the  scattering  volume, 
Fig.  11a  and  b,  appears  completely  random  relative  to  the  prep  shots. 

This  might  be  expected  since  the  beam  axis  moved  around  from  shot  to  shot 
relative  to  that  volume.  However,  the  beam  axis  was  always  above  the 
lower  half  of  the  scattering  volume,  and  there  beam  shot  line  centers 
were  very  consistent,  as  shown  in  Figs.  11c  and  d.  (Note  that  the  shift 
in  average  line  center  for  prep  shots  included  in  Fig.  11  relative  to  the 
full  scattering  volume  prep  shots  in  Fig.  10  is  believed  to  be  due  to 
differences  in  the  upper  and  lower  halves  of  the  fibre  bundles.  Thus, 
it  is  the  relative  shift  between  beam  and  prep  shots  which  is  of  interest 
here . ) 
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We  note  that  a 3 A line  center  shift  implies  a drift  velocity 
of  1.3  X 10  cm/sec.  At  10  cm  density,  this  would  imply  a plasma  elec- 
tron current  of  > 2 kA/cm  . This  is  an  order  of  magnitude  higher  than 
we  would  expect  from  plasma  diamagnetic  currents,  or  residual  net  cur- 
rent density  after  beam  passage  through  the  plasma.2'  However,  any  lack 
of  symmetry  in  the  interaction,  or  residual  radial  electric  fields, 
might  cause  gross  plasma  motion  which  would  not  have  currents  associated 
with  it.  Extreme  asymmetries  of  the  type  suggested  by  VanDevender,  et 
al.,  1 as  an  explanation  for  apparent  drift  velocities  of  ~ 10“  cm/sec 
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In  this  equation  is  the  sum  of  electron  and  ion  transverse  energies 

per  particle,  Bq  is  the  applied  magnetic  field,  V,^  is  the  observed 

voltage,  T is  the  diamagnetic  loop  circuit  integration  time  constant, 

and  R^  and  are  the  loop  and  conducting  wall  radii,  respectively.  The 

average  signal  V.^  was  9 V (with  a 20$  standard  deviation),  t = 5. 6 us, 

B = 0.4  W/m2  and  (l  - R?/R2)  = 0.4.  Taking  the  beam  cross  sectional 
o L w 

area  of  40  cm2  gives  A(n  e^)  = 6 x 101'  eV/cm3,  or  60  eV  per  electron 
ion  pair  at  1015/cm3  density.  However,  we  will  very  shortly  see  that  at 
least  100  cm2  is  a more  reasonable  estimate  of  the  heated  cross-section- 
al area,  giving  24  eV  per  electron  ion  pair  at  101  ~7cm° . From  Figs. 

5 and  8 we  see  that  the  electron  temperature  rise  according  to  Thomson 
scattering  is  only  ~ 5 eV.  Even  assuming  equal  energies  in  electrons 
and  ions,  the  discrepancy  here  is  a factor  of  4,  remarkably  close  to  the 
factor  of  6 observed  by  Arzhannikov,  et  al . ,x  at  ~ lO^/cm  density. 

This  may  be  due  to  energy  in  rotational  drift  motion  which  was  noted  as  a 
possible  explanation  for  the  larger  than  expected  liner  center  shifts. 

The  presence  of  an  energetic  tail  on  the  electron  distribution  function, 
as  predicted  by  e-e  instability  theory,  0,3  may  also  explain  the  dis- 
crepancy. Finally,  it  may  be  the  residual  effect  of  beam  diamagnetism.0' 
If  it  is  a tail,  it  must  be  of  order  500  eV  or  more  in  order  not  to 
thermalize  with  the  main  distribution  in  500  - 400  ns  by  collisions. 

The  heated  plasma  area  was  obtained  in  a separate  experimental 
run  in  which  a 40  - 50  kA  peak  current,  900  keV  peak  voltage,  60  ns  beam 
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was  injected  into  the  theta  pinch.  The  plasma  density  was  (l  ± .5)  X 
1015/  cm  and  the  temperature  ~2.5  eV.  __  The  laser  scattering  system 
was  not  in  operation  during  this  run.  These  values  for  density  and 
temperature  were  obtained  in  a series  of  prep  shots  taken  after  the  probe 
shots  were  completed.  Prep  shots  were  taken  both  with  and  without  the 
probe  housing  (and  probe)  in  place.  There  was  no  significant  difference.^ 
The  magnetic  probe,  as  described  in  Sec.  Ill,  located  m into  the 
theta  pinch  (and  1/2  m in  front  of  the  second  diamagnetic  loop)  was  used 
to  measure • local  magnetic  field  changes  resulting  from  the  beam  plasma 
interaction.  Figure  12a  shows  the  interaction  geometry  at  the  probe 
port  for  this  run.  (The  beam  position  was  obtained  using  solid  targets 
in  the  beam  path  at  the  probe  port  as  well  as  by  x-ray  pinhole  photo- 
graphy on  the  probe  housing  and  on  a coarse  grid  of  Tungsten  wires 
placed  at  the  same  axial  position.)  Figure  12b  shows  a typical  Faraday 
cup  waveform  for  time  reference.  The  results  for  the  change  in  axial 
magnetic  field  (ABz)  and  the  horizontal  field  (B^)  are  shown  in  Figs. 

13a  and  13b.  Note  that  in  an  axisymmetric  experiment,  B^  would  be  the 
azimuthal  field  component,  Bp . B^  was  obtained  from  the  probe  signal, 
dB^/dt,  by  the  standard  technique  of  a passive  RC  integrator  (having  a 
2 >;s  time  constant)  at  the  osci lloscope . However,  it  was  necessary  to 
display  dB^/dt  and  graphically  integrate  the  oscillograph  to  obtain  ABz 
because  of  the  voltage  induced  in  the  probe  by  the  theta  pinch.  These 
curves  have  been  published  in  our  previous  article, but  our  discussion 
here  will  be  greatly  facilitated  by  having  them  at  hand. 

At  the  moment  our  interest  is  the  diamagnetic  area  in  Fig.  13a 
after  the  beam  has  passed.  (A  diamagnetic  ABz  signal  is  negative  in  the 
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figure.)  We  note  that  it  is  considerably  larger  than  the  region  in  which 
current  is  flowing,  as  indicated  by  the  signals  in  Fig.  ljb.  More- 
over, since  the  beam  geometry,  obtained  from  targets  and  x-ray  pinhole 
photography  (Fig.  12a),  and  the  signals  both  indicate  that  the  probe 
is  measuring  ABz  across  a chord  in  the  heated  plasma  rather  than  a 
diameter,  we  infer  a disturbed  area  of  100  - 125  cm2  after  the  beam 
pulse  (see  the  profiles  at  125  or  175  ns). 

The  100  gauss  depth  of  the  diamagnetic  well  implies  a A(  e ) of 
~ 20  eV  per  electron-ion  pair  for  1015/cm"  density.  Diamagnetic  loop 
amplitudes  corresponding  to  about  half  this  temperature  change  were  ob- 
served on  these  shots  (assuming  100  cm^  heated  cross  section).  These 
diagnostics  imply  ~ 50  j/m  of  beam  energy  deposited  in  the  plasma,  assum- 
ing isotropy,  for  a coupling  efficiency  of  5^/m. 

It  is  interesting  to  note  that  the  disturbed  cross-sectional 
area  indicated  in  Fig.  13a  even  as  early  as  50  ns  is  considerably  larger 
than  the  area  in  which  the  beam  current  (and,  therefore,  the  plasma  cur- 
rent) is  flowing.  The  speed  of  cross  field  energy  transport  implied  by 
this  is  ^ 40  cm/us,  a value  considerably  larger  than  that  obtained  in  an 
ordinary  turbulent  heating  experiment  by  Aranchuk,  et  al . Thus,  not 
only  does  beam  heating  of  a plasma  avoid  skin  effect  difficulties  of 
ordinary  turbulent  heating,  it  also  rapidly  heats  the  surrounding  volume, 
possibly  by  wave  energy  transport.  78 

Turning  now  to  the  profiles,  there  are  several  character- 
istics of  interest.  Firstly,  as  already  noted,  the  net  current  implied 
by  this  profile  during  the  beam  pulse  locates  the  beam  in  a position  con- 
sistent with  target  and  x-ray  diagnostics  (i.e.  above  the  discharge  tube 
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axis  by  ~ 2 cm).  Secondly,  the  net  current  position  shifts  from  2 - 3 cm 
off  axis  to  very  nearly  on  axis  at  the  end  of  the  beam  pulse.  The  reason 
for  this  is  not  known,  although  it  may  be  reflecting  a movement  of  the 
beam  as  its  current  drops  from  its  maximum  to  zero  at  the  end  of  the 
pulse.  The  magnitude  of  the  net  axial  current  is  also  of  interest. 

Taking  into  account  the  system  geometry  (Fig.  12a),  ~ 500  A of  net  cur- 
rent is  implied,  i.e.,  1^/100.  This  is  a factor  of  3 greater  than  would 
be  predicted  by  sharp  beam  boundary  theory,27’" 1 and  a factor  of  1000 
greater  than  for  the  more  realistic  beam  radial  density  profiles  of 
Kuppers,  et  al.,40  ignoring  return  current  damping."7  Assuming  classi- 
cal resistivity,  the  damping  time  T (see  Sec.  IIA)  for  a 3*5  cm  radius, 
2.5  eV  plasma  is  7 Lis.  Thus,  a net  current  density  of  ~ 1$  might  be 
expected  in  ~ 70  ns.  The  fact  that  it  appears  virtually  instantaneously 
(by  the  end  of  the  beam  rise  time),  and  then  changes  very  slowly  during 
the  main  part  of  the  beam  pulse  appears  to  be  consistent  with  a collision 
frequency  > 10  times  classical  at  first,  and  perhaps  2-3  times  classi- 
cal during  the  next  50  ns.  (T  increases  to  ~ 20  as  at  5 eV.)  These 
collision  frequencies  are  consistent  with  values  obtained  by  calculating 
the  heating  rates  implied  by  the  magnetic  profiles." 


B . Partially  Ionized  Case 

The  use  of  a partially  ionized  plasma  adds  the  possibility  of 


substantial  ionization  and  radiative  energy  losses  to  the  other  processes 
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helium.  In  fact  , twice  these  minimum  energies  are  required  since  line 
radiation  energy  losses  must  be  considered.  In  previously  reported  re- 
sults^1 from  the  present  experiment,  a ~ 30  kA,  550  keV,  70  ns  beam  was 
injected  into  a 5 X lO^/cm3  density  helium  plasma  with  neutral  helium 
density  n^  S' 3-5  X lO^/cnr  (100  mTorr  pressure).  The  applied  magnetic 
field  throughout  the  experimental  system  was  5 kG  in  this  run,  as  well  as 
in  the  partially  ionized  cases  to  be  discussed  below.  The  plasma  density 
and  temperature  observed  after  beam  passage  were  ~ 7 X lO^/cm3  and  up  to 
5 eV.  Using  these  numbers  and  40  eV  to  produce  each  electron- ion  pair  to 
estimate  the  energy  input,  we  obtain  n (5/2  kTg  + 4o)/70  ns  « 5 X 1023 
eV/cnT  - sec  as  the  average  energy  deposition  rate.  This  represents  a 
coupling  efficiency  of  ~ J^/m  assuming  a heated  plasma  cross  section 
equal  to  the  beam  area  (20  cm2).  If  the  heated  area  is  2 V 2 that  area 
as  was  the  case  for  the  highly  ionized  experiments  discussed  in  Sec.  IIA, 
the  implied  coupling  efficiency  is  ~ 7$/m.  (The  factor  of  2 1/2  is  only 
conjectural  in  the  partially  ionized  experiments  since  no  magnetic  probe 
scans  were  taken.) 

Two  additional  experimental  runs  were  taken  in  which  the  beam 
was  injected  into  partially  ionized  helium  at  100  mTorr  pressure.  In 
the  first,  a 1 MeV,  50  kA  peak  current,  5 cm  diameter  beam  of  70  ns 
duration  was  injected  into  an  8 x 10l3/cm3  density,  2 eV  plasma  (the 
"8  x lo13  case").  In  the  second  run,  the  beam  was  900  keV  and  40  kA  peak 
current,  and  the  initial  plasma  density  was  4 x 1014/cm°,  with  the  re- 
maining parameters  being  the  same  (the  "4  x 1014  case").  Figures  14  and 
15  present  the  plasma  electron  temperature  and  density  obtained  by 
Thomson  scattering.  The  time  interval  covered  was  70  - 450  ns  after 
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beam  front  arrival  at  the  scattering  port,  as  defined  in  Fig.  ba.  (As 
discussed  in  Sec.  Ill,  scattering  during  the  beam  pulse  was  not  possible 
due  to  the  high  level  continuum  light  during  the  beam  pulse.) 

The  data  in  Figs.  1^-  and  15  both  show  the  following  character- 
istics: 1.)  The  plasma  density  rises  from  its  initial  value  to  ~ l015/cm3 

by  the  end  of  the  beam  pulse,  and  2.)  the  temperature  is  6 - 8 eV  just 
after  the  beam,  and  falls  over  ~ 100  ns  to  about  3 eV  where  it  stabilizes. 

These  two  runs  give  average  energy  deposition  rates, 
n^(3/2  kTe  + 4o)/70  ns,  slightly  higher  than  the  previously  reported 
case,24  namely  (7  - 12)  X 102-  eV/cmJ  - sec.  Because  the  injected  beam 
energy  for  these  runs  was  ~ 3 times  that  of  the  previous  case,  the  re- 
sulting coupling  efficiencies  for  these  higher  density  cases  were  lower, 
namely  ~ 2/5/m  assuming  only  the  beam  area  is  heated,  and  ~ 5$/m  if  2 1/2 
times  that  area  is  heated,  as  previously  discussed. 

In  order  to  understand  these  results  a one  dimensional  ioniza- 
tion and  heating  model,  similar  to  that  used  previously24  was  constructed. 

It  is  discussed  in  detail  in  Appendix  B.  To  summarize,  energy  is  de- 
posited resistively  in  the  plasma  at  a rate  T|j'  , where  T|  is  the  plasma 
resistivity  and  j is  the  return  current  density  in  the  plasma  (assumed 
equal  to  the  beam  current  density).  The  resulting  changes  in  plasma 
density  and  temperature  are  followed  in  time  by  solving  a coupled  set  of 
differential  equations  Eqs.  ^Bl  - B7))  for  the  densities  of  neutral, 
singly,  and  doubly  ionized  helium,  and  for  the  temperature  of  these 
species.  The  principal  energy  loss  processes  included  are  ionization 
and  line  radiation,  although  Eq . (B5)  for  the  electron  energy  contains 
several  others.  Rate  coefficients  are  temperature  dependent  as 
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appropriate ' c ’ 1 (Eqs.  (Bll)  and  (B12)).  The  quantity  T|  includes  classi- 
cal (Eq.  B8),  as  well  as  anomalous  resistivities  due  to  ion  acoustic  and 
e-e  mode  turbulence  (Eqs.  B9  and  BIO),  which  are  discussed  in  Appendix  B, 


The  model  ignores  such  plasma  dynamic  and  kinetic  effects  as 
expansion  and  end  loss.  This  is  equivalent  to  assuming  the  plasma  to  be 
spatially  uniform.  This  should  be  consistent  with  the  low  observed 
temperatures,  the  short  timescale  of  interest,  and  classical  transport 
and  thermal  expansion. 

The  plasma  current  density  is  assumed  to  have  a 10  ns  e- fold 
rise  time,  a 70  ns  width  ( fwhm) , and  a 10  ns  e-fold  fall  time. 

Results  obtained  from  this  model  are  shown  together  with  the 
experimental  data  in  Figs.  14  and  15.  We  see  that  the  theoretical  plasma 
temperature  shoots  up  at  early  time  (Figs.  14a  and  15a).  This  is  because 
there  are  relatively  few  electrons  to  share  the  energy  input.  (This  is 
enhanced  by  the  inverse  dependence  of  the  anomalous  resistivities  on 
density,  Eqs.  (B9)  and  (BIO).)  The  exponential  dependence  of  ionization 
on  temperature,  Eq . (Bll),  means  that  only  when  the  temperature  is  above 
~ 10  eV  will  rapid  ionization  occur.  This  occurs  in  less  than  10  ns  in 
the  8 x 1013  case,  and  the  density  rises  rapidly,  doubling  in  25  ns  and 
redoubling  in  less  than  50  ns  in  this  case  (Fig.  l4a ) . The  rate  of 
energy  input  then  decreases  (since  the  anomalous  resistivities  decrease) 
and  the  energy  is  divided  among  more  electrons,  depressing  the  tempera- 
ture. By  the  end  of  the  beam  pulse  in  the  8 x lO1^  case,  the  temperature 
is  only  V5  of  its  maximum  value.  At  this  time,  the  energy  input  stops 
(since  j decreases  to  zero),  the  temperature  rapidly  decreases  to  less 
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than  10  eV  and  further  ionization  ceases.  In  the  4 x 1014  case,  the 
temperature  peaks  at  a much  lower  value.  Therefore,  the  ionization  pro- 
ceeds at  a slower  rate  than  in  the  first  case,  and  the  density  when  the 
beam  (and  heating)  pulse  is  over  is  only  three  times  the  initial  value. 

We  note  that  the  plasma  conditions  after  the  beam  pulse  are  insensitive 
to  small  changes  in  the  initial  conditions.  For  example,  decreasing  the 
initial  electron  density  to  5 X 1014/cm3  changed  the  density  by  2$  and 
decreased  the  temperature  by  Vp  at  t = 100  ns.  This  was  in  spite  of  the 
fact  that  the  maximum  temperature,  at  t = 35  - 40  ns,  increased  by  15^-. 

Comparing  the  theoretical  and  experimental  results,  we  see 
that  the  theoretical  density  rises  too  quickly  in  the  8 x lO1^  case  and 
too  slowly  in  the  4 x 101"  case.  Agreement  with  the  final  value  is  quite 
good  in  the  lower  density  case,  but  it  is  low  in  the  higher  density  case. 
In  both  cases,  the  experimental  temperatures  appear  to  drop  to  their 
asymptotic  values  (~  3 eV  in  both  cases)  much  more  quickly  than  the  model 
predicts.  This  is  probably  a result  of  either  nonclassical  energy  trans- 
fer to  ions  (possibly  due  to  the  presence  of  ion  acoustic  turbulence),  or 
nonclassical  radial  energy  transport  (as  by  waves).  For  example,  if  half 
the  plasma  energy  at  ~ 100  ns  is  apportioned  to  ions  (which  are  ^ 1 eV 
from  classical  heating),  temperature  agreement  would  be  quite  good  at  that 
time.  Other  processes  for  plasma  energy  loss  (such  as  via  impurity 
radiation,  assuming  of  order  1$  Nitrogen,  Carbon  and/or  Oxygen)  are  too 
slow  to  account  for  the  necessary  electron  cooling  rate. 

Since  we  have  used  only  resistive  heating  in  our  model,  it  is 
reasonable  to  ask  what  will  happen  to  these  results  if  nonresistive 
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heating  by  the  electron-electron  instability  is  added.  In  Sec.  V,  we 
shall  see  that  direct  heating  from  this  instability  can  be  expected  to 
be  greatest  near  the  beginning  of  the  beam-plasma  interaction  for  our  ex- 
perimental conditions.  Therefore,  in  the  4 x 1014  case,  we  have  arbi- 
trarily increased  Q in  Eq.  (B5)  by  1024  eV/cm3  - second  for  10  ns  from 
t =6  to  16  ns.  The  resultant  density  at  100  ns  was  10$  higher,  and  the 
temperature,  although  higher  at  25  ns,  was  virtually  the  same  (2 $ lower) 
at  100  ns.  We  have  not  carried  the  numerical  calculation  beyond  200  ns 
since  its  1-dimensionality  and  lack  of  plasma  transport  mechanisms  limit 
its  validity  to  short  times. 

A final  numerical  result  to  be  noted  is  the  deposited  energy 
per  electron- ion  pair.  From  Fig.  14  ( 8 x 1013  case),  the  plasma  energy 

density,  5/2  n kT  , is  ~ 1.5  X 101  eV/cm3.  An  energy  of  5-^  X 101 

P ® 

eV/cnf  was  deposited  ( resistively ) in  that  numerical  run.  Thus,  approxi- 
mately 4 x 101  eV/cm  remains,  implying  40  eV  was  required  to  produce 
each  electron- ion  pair  , as  previously  assumed.  (The  approximate 
equality  of  the  ionization  and  excitation  rate  coefficients,  Eqs.  (B10^ 
and  (Bl2),  explains  the  need  for  about  1.6  times  the  24.  eV  ionization 
energy  for  helium.) 

In  addition  to  Thomson  scattering  measurements,  we  have  ob- 
served microwave  emission  in  X-band  (8-12  GHz)  and  optical  emission 
of  helium  I and  II  lines,  and  have  made  diamagnetic  loop  measurements. 

The  microwave  measurements  revealed  strong  emission  near  the  relativistic 
cyclotron  frequency  corresponding  to  the  diode  voltage  and  the  applied 
magnetic  field  (5  kG).  This  emission  occurred  at  the  beginning  of  the 
beam  pulse  and  had  K 20  ns  full  width  at  half  maximum.  The  optical  line 
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emission,  however,  could  not  be  accounted  for  by  assuming  all  the  ener^ 
deposition  occurred  during  this  time.  A more  uniform  deposition  rate 
was  required.'  This  leads  to  the  conclusion  that  the  X-band  radiation 
was  not  associated  with  the  bulk  of  the  energy  deposition. 

Taking  up  the  diamagnetic  loop  results  now,  Fig.  1 6 shows  a 
typical  diamagnetic  loop  signal  (after  subtracting  out  the  uncompensated 
portion  of  the  induced  voltage  from  the  theta  pinch)  for  the  X 101 
case,  (it  is  representative  of  the  other  case  as  well.)  Very  similar 
amplitude  and  shape  signals  were  obtained  on  the  two  diamagnetic  loops 
when  both  were  used,  indicating  the  uniformity  of  the  interaction  over 
the  4 m length  of  the  interaction  region  for  this,  the  partially  ionize< 
case,  as  it  was  in  the  highly  ionized  case.  (The  only  difference  was 
that  the  diamagnetic  loop  signal  at  the  upstream  end  of  the  system  some- 
times showed  oscillations  of  the  type  to  be  discussed  below.’'  Because 
of  the  slow  response  time  of  the  diamagnetic  loop  circuit  the  signal 
voltage  is  not  simply  related  to  the  instantaneous  plasma  energy  density 
as  measured  by  Thomson  scattering.  We  assume  the  theoretical  density  a: 
temperature  time  histories  as  shown  in  Fig.  1‘-  to  be  correct,  and  appro: 
imate  the  product  n^kT^  by  a triangular  shaped  transverse  energy  pulse 
of  full  width  100  ns  and  peak  of  5 x 101  eV/cm  at  t = 50  ns.  We  then 
calculate  the  expected  oscilloscope  voltage  as  a function  of  time  V 
using  the  circuit  shown  in  Fig.  17a  to  approximate  the  diamagnetic  loop 
circuit.  With  L = 22  uh  and  RC  = 10  as,  the  result  is  given  in  Fig.  171 
The  signal  peak  is  only  a factor  of  two  smaller  than  the  observed  signa 
Fig.  1 . 
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In  Fig.  17,  we  calculated  the  expected  diamagnetic  loop  signal 
from  a triangular  shaped  (in  time)  pulse  of  diamagnetism  which  was  chosen 
as  an  approximation  to  the  plasma  diamagnetism  implied  by  Fig.  14.  How- 
ever, it  could  just  as  well  represent  beam  diamagnetism  since  n^  « 4 x 
10l;L/cm  (2  kA/cm2)  would  require  only  150  keV  transverse  energy  per  beam 
electron  to  produce  the  6 y 101  eV/cm  peak  we  used  in  Fig.  1J.  In  the 
final  series  of  shots  to  be  discussed  in  this  paper,  12  cm"  and  40  cm" 
beams  at  ~ 900  keV  were  injected  into  partially  ionized  hydrogen  at 
^ 5 X 101 J /cm  and  ~ 1 eV.  The  beam  current  density  was  in  the  range 
~ 1.2  to  ~ 5*5  kA/cm"  for  the  smaller  area  beam  and  was  ~ 1.2  kA/cm"  for 
the  larger  beam.  (Density  and  temperature  as  a function  of  time  from 
laser  scattering  gave  results  quantitatively  similar  to  the  partially 
ionized  helium  data  shown  in  Figs.  14  and  15.  There  were,  however,  too 
few  shots  at  any  given  beam  condition  to  draw  a graph  similar  to  those 
figures  for  these  shots.)  Figure  1 shows  typical  diamagnetic  signals 
for  these  shots.  In  Fig.  19a  we  plot  the  amplitude  of  the  diamagnetic 
loop  near  the  calorimeter  against  the  calorimeter  energy.  The  smaller 
area  beam  shows  a clear  trend  toward  loop  signal  being  proportional  to 
propagated  beam  energy.  (Note:  since  the  beam  voltage  pulse  duration 

was  virtually  the  same  for  all  shots  plotted,  only  the  current  varies 
to  produce  the  variation  in  beam  energy.)  The  larger  area  beam  does  not 
show  this  trend. 

In  Fig.  19b,  tipkT  obtained  from  laser  scattering  is  plotted 
against  diamagnetic  loop  amplitude  for  those  shots  from  Fig.  19a  for 
which  scattering  data  was  taken.  There  is  some  indication  here  that 
plasma  energy  after  the  beam  pulse  and  the  diamagnetic  loop  signal  are 
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correlated.  As  in  the  helium  cases,  the  diamagnetic  loop  is  indicating 
more  than  an  order  of  magnitude  more  transverse  energy  than  the  laser 
scattering,  presumably  for  the  same  reasons.  The  data  of  Fig.  19a 
suggests  that  beam  diamagnetism  cannot  be  ruled  out  as  a contributing 
factor,  as  least  for  the  smaller  area  beam.  However,  an  equally  allowab’ 
explanation  is  stronger  beam-plasma  interaction  for  the  higher  current 
density  beams. 

Returning  to  Fig.  18,  we  note  that  both  signals  show  oscilla- 
tions with  the  same  period,  probably  the  magnetosonic  oscillations  pre- 
viously discussed. : ,;l"  »44»  ' The  more  prominent  oscillations  on  th< 

diamagnetic  loop  signals  here  as  compared  to  the  helium  shots  (e.g.  - 
Fig.  1 6)  may  be  a result  of  the  magnetosonic  waves  being  damped  by  more 
collisions  in  the  helium  case/1'"  In  addition  we  see  that  the  amplitudes 
are  virtually  the  same,  again  indicating  the  uniformity  of  the  inter- 
action of  the  entire  system  length. 

V.  THEORETICAL  CONSIDERATIONS 

In  this  section  we  address  the  questions  of  what  beam- to-p lasma 
energy  transfer  mechanisms  are  expected  under  the  conditions  of  our  ex- 
periments, and  how  much  energy  we  expect  the  beam  to  lose  according  to 
the  appropriate  theory.  We  first  present  a qualitative  discussion  of  tl 
most  probable  loss  mechanisms  (direct  electron-electron  two  stream  and 
return  current  heating),  including  their  time  dependent  characteristics 
These  are  subsequently  applied  to  our  experiments. 

According  to  linear  theory  two  types  of  waves  are  expected  durinj 
IREB  - plasma  interaction.  The  first  type  corresponds  to  Langmuir  wave; 

kk 


w ' 


excited  by  the  electron-electron  two  stream  instability  (e-e  mode)  with 


wavelength  of  the  order  c/u^.  We  refer  to  the  energy  density  of  these 


waves  as  W . The  second  type  corresponds  to  ion  acoustic  waves,  W^, 


which  can  be  excited  when  the  induced  return  current  drives  electron- ion 
instabilities  (e-i  modes).  An  approximate  criterion  for  the  excitation 


of  such  instabilities  is  v,/c  » T./T  (v,  > c if  T.  « T ),  where 

ds  led  s 1 e 


c = 
s 


k(T_  + T^)M  j ^ and  vA  is  the  plasma  electron  drift  velocity  rela- 


e 1 J d 

tive  to  the  mass  M ions.  Two  energy  transfer  mechanisms  are,  therefore. 


possible.  One  is  the  direct  interaction  of  the  beam  electrons  in  reso- 
nance with  the  waves  W (i.e.,  x - c (k  »j>/|jg>J)  = 0 where  £ is  the  b 


earn 


electron  momentum  vector).  The  second  is  scattering  of  the  plasma  elec- 


trons forming  the  return  current  on  W , which  results  in  anomalous lv  high 

a 


resistance.  In  order  to  compute  the  energy  transfer  rates,  the  wave 


energy  levels  W and  W must  be  known.  The  computation  of  these  levels 

1 a 


has  been  the  most  controversial  aspect  of  IREB  plasma  heating. 

An  extensive  amount  of  work  has  focused  on  the  determination  of 


W and  the  associated  energy  coupling  length  (£)  on  the  basis  of  con- 


or.  po  o T 

vective  quasilinear  theory.  9 9 Such  considerations  produced  the 


length  given  by  Eq.  (4).  As  discussed  in  Sec.  II  experiments  have  shown 
this  to  be  much  too  short.  This  is  physically  expected  since  for  any 


parameters  of  interest,  the  magnitude  of  W violates  the  assumptions  of 


quasilinear  theory.  An  alternative  possibility  is  that  the  amplitude 


W is  limited  by  nonlinear  wave-wave  interactions.  The  general  idea  of 


this  concept  is  that  beyond  a certain  level  of  W , wave  energy  is  trans- 


ferred into  a nonresonant  region.  A stationary  state  can  then  be  achieved 


j 


^5 


where  the  energy  transfer  from  the  beam  to  the  (resonant)  Langmuir  waves 
is  balanced  by  the  transfer  of  wave  energy  to  the  nonresonant  region. 

As  a result,  W is  maintained  at  a sufficiently  low  level  while  dissipa- 
tion occurs  only  in  the  nonresonant  regions.  The  energy  loss  rate  for 
the  beam  electrons  is  then  proportional  to  2 6 W , where  6 is  the  insta- 
bility growth  rate  given  by  Eq . (3)-  Early  attempts  to  apply  these  con- 
cepts to  IREB  plasma  experiments'^  ’ 1 were  confined  within  the  frame- 
work of  validity  of  weak  turbulence  theory  (i.e.  the  real  part  of  the 
frequency  obeys  the  linear  dispersion  relation).  These  efforts  still 
failed  to  reconcile  the  differences  between  theory  and  experiment.  How- 
ever, Papadopoulos  3 showed  that  for  any  reasonable  parameters  applicable 
to  present  day  intense  beam  plasma  interaction  experiments  the  weak 
turbulence  theory  is  not  valid  and  inclusion  of  nonlinear  frequency 

shifts  is  important.  It  was  shown  that  when  W /n  kT  > (v  /c)2,  the  wave 

1 p e e 

spectrum  becomes  unstable  to  a secondary  instability  similar  to  the 
oscillating  two  stream  instability  (OTSl)?  ,72  (Within  this  context  the 
instability  is  known  also  as  the  modulational  or  modified  decay  insta- 
bility.) This  process  transfers  energy  to  electron  plasma  waves  with 
lower  phase  velocities  (shorter  wavelengths)  and  associated  low  frequency 
ion  waves.  The  lower  phase  velocity  plasma  waves  can  be  linearly 
Landau  damped  by  the  tails  of  the  plasma  electron  distribution  function. 
Note  that  these  wave  processes  viewed  in  configuration  space  correspond 
to  plasma  waves  trapped  in  low  density  regions  and  have  been  given  the 
name  of  plasma  solitons,  cavitons  and  spikons.  It  has  been  shown74 
that  they  are  equivalent  representations  of  the  OTSI.  It  should  also  be 
noted  that  among  the  new  concepts  introduced  by  the  strong  turbulence 
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theory  is  the  possibility  of  exciting  ion  waves  W created  directly  by 

SL 

the  ponderomotive  force  exerted  on  the  plasma  by  W even  in  the  absence 

of  return  current  driven  instabilities  ( e- i modes). 

On  the  basis  of  this  model,  the  details  of  which  are  available 

elsewhere,  the  time  sequence  of  events  is  as  follows.  Upon  injection 

of  the  beam  in  the  plasma  the  waves  W grow  rapidly  until  they  reach  a 

1 

level  such  that  their  removal  from  the  instability  region  is  faster  than 
the  e-e  instability  growth.  This  is  given  by  (see  Eqs.  (25-2o)  of  Ref. 

53) 
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where  M is  the  ion  mass.  They  subsequently  decay  to  Langmuir  waves  non- 
resonant with  the  beam  electrons  (w  ) and  modified  ion  plasma  waves  (w  ). 

2 a 

This  is  illustrated  in  Fig.  20.  A quasistationary  state  can  be  estab- 
lished on  the  basis  of  the  following  arguments:  The  presence  of  finite 
amplitude  ion  waves  Wa  can  produce  an  ac  resistivity  for  waves  near  the 
plasma  frequency  with  effective  collision  rate  as  discussed  by  Dawson 

and  Oberman7  and  a dc  resistivity3  ' ' with  rate  vw  = k>  v*,  where 

J o D H 

kip  is  the  characteristic  wave  number  Debye  length  product  for  the  ion 
waves.  When 


VH  " 2 


(iM 


k7 


the  e-e  instability  is  stabilized  and  the  energy  deposition  rate  via 


waves  de  /dt  will  be  given  by 
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de 


1 

dt 
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where  D is  the  dielectric  function  of  the  plasma.  In  addition  to  this, 


the  presence  of  the  dc  collision  frequency  will  provide  an  energy 


deposition  mechanism  due  to  the  return  current  j = n^ev^.  This  will  be 
given  by 


de 


2 _ -,#.2  = 


it-TTv* 


J = 


kw 
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dt 
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Notice  that  if  the  Coulomb  collision  rate  v . is  greater  than  v!t  or  v*  , 

cl  Ho 

then  it  replaces  them  in  these  arguments.  Furthermore  if  an  e- i insta- 
bility is  present,  W and  therefore  the  effective  collision  frequencies 

cl 

o Q 

might  have  to  be  determined  by  other  considerations.  The  energy  levels 

W and  W are  determined  by  the  condition  that  the  ion  waves  are  margin- 

1 2 

ally  stable.  This  gives  L_Ref.  35  Eqs.  (39-^MH 


n kT  0)  u) 

p e p p 
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It  should  be  noted  that  these  relations  have  been  verified  by  computer 
simulations  using  particle  and  mode  coupling  codes.7. 
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We  now  proceed  to  apply  these  concepts  to  the  present  experiment. 

In  order  to  be  more  precise  quantitatively  we  select  first  the  case 
where  the  IREB  ( n^  « 5 X 1011/cm  ) was  injected  into  a highly  ionized 
plasma  with  density  (7  ± 3)  X IQ1 4 /cm"'  and  initial  temperature  « 

3 eV . For  these  parameters,  e- i instabilities  are  not  expected,  since 
v^  « v . We  consider  the  beam  pulse  shape  shown  in  Fig.  21  and  see  if 
Eqs.  (13- 17)  predict  the  energy  absorption  data  and  the  field  penetration 
timescale  as  measured  by  the  laser  scattering  and  magnetic  probes.  The 
time  t = 0 is  the  time  the  beam  front  arrives  at  the  particular  diag- 
nostic port.  We  compare  the  observation  first  with  the  approximate 
analytic  results,  and  then  the  results  of  a numerical  solution  of  the 
nonlinear  equat ions . r7 

For  the  parameters  of  this  experiment  and  for  the  energy  deposi- 
tion during  the  initial  stage  we  find  from  Eq . ( 13T~th'at  W max/n  kT  = 

l p e 

TTlcl  X / 

0.32.  The  actual  energy  loss  of  the  beam  is  given  by  Ae  = W ( 1 + 

o 1 

« 2 W maX  « I.36  x 1015  eV/cm:.  This  energy  is  delivered 
1 

■s 

during  a 5 ns  time  interval  around  t = 20  ns,  since  the  instability  does 

not  start  (2  6 < v „)  until  t ss  15  ns.  At  this  time  both  v*  and  v*  reach 
c/j  Ho 

a maximum  and  then  relax  towards  their  quasistationary  values.  For  the 

case  under  consideration  these  are  given  by  W /n  kT  « W /n  kT  « 8 6/u)  « 

1 P e 2Pe  p 

2 x 10“r , Wa/npkTg  * 10'4  and  kX^  w .1  - .2  ( for  y = 3 and  5 % 1/3 ) . 

With  the  above  values  we  find  that  for  the  rest  of  the  beam  pulse  the 
energy  deposition  rate  due  to  wave  damping,  given  by  Eq.  (15)»  is 


dxD(u),  k) 
3uu 


while  Che  one  due  to  the  return  current,  Eq . (lo,  is 


de 

2 

dt 
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x 1023 


eV 


cm  - sec 


The  total  energy  lost  by  the  beam  according  to  this  model  is  ~ 

3.5  X 101'"  eV/cm‘  . For  a 40  cm^  beam,  this  gives  ~ 1.5  X 1016  eY/cm 
during  the  entire  pulse.  This  is  to  be  compared  with  ~ 5 X lO17  eV/cm 
deposited  energy  measured  by  laser  scattering  assuming  100  cm2  heated 
plasma  cross  section,  (Fig.  6a)  and  ~2  x 1016  from  the  diamagnetic  loop 
average . 

We  have  checked  the  conclusions  of  the  simplified  analytic  model 
by  numerically  solving' ' the  exact  mode  coupling  equations  which  are 
derived  in  Ref.  (33)>  including  the  effect  of  finite  beam  rise  time  and 
classical  collisional  damping.  The  results  for  the  energy  deposition, 
low  frequency  resistivity,  and  wave  spectrum  as  a function  of  time  are 
shown  in  Figs.  22-2<- . The  analytic  results  are  included  for  comparison 
and  are  seen  to  be  consistent  with  the  discussion  we  have  presented.  Note 
that  the  computational  results  for  energy  deposition  (Fig.  22)  do  not  in- 
clude resistive  heating  contributions.  From  Fig.  23  we  can  see  an  early 
time  resistivity  which  is  more  than  an  order  of  magnitude  larger  than 
classical,  and  a later  time  resistivity  of  about  twice  the  classical 
value,  as  required  to  explain  the  magnetic  probe  results  (Fig.  13).  The 
high  and  low  frequency  wave  spectra  during  the  various  stages  of  the 
interaction,  shown  in  Fig.  24,  demonstrate  that  modes  with  phase  velocity 
in  resonance  with  the  beam  (marked  by  +)  dominate  only  at  early  time. 

In  the  higher  density  highly  ionized  experiments  for  which  re- 
sults are  shown  in  Fig.  6 Zi-e.,  n = (1.5  i .5)  X 1015/cma  and  n = 

P p 

(3.1  i 1)  X 10"_/crr^,  2 5 < , for  all  time  during  the  pulse  assuming 
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that  T = 5 eV  and  9 as  1/3.  Therefore,  we  do  not  expect  the  e-e  in- 
e 

stability  to  be  excited.  For  the  case  (3.1  ± 1)  X lO1^  cm',  as  shown  in 
Fig.  8c,  classical  heating  by  dissipation  of  the  return  current  can 
account  for  the  observed  increase  in  temperature.  However,  this  is  not 
true  for  the  (1.5  ± *5)  X id1"/  cm3  case.  At  this  stage  we  can  only 
speculate  that  perhaps  some  nonuniformity  in  the  beam  current  density 
can  produce  a sufficient  local  heating  at  early  time,  to  allow  2 6 > 

and  the  e-e  instability  to  be  excited. 

We  proceed  now  to  examine  the  partially  ionized  results.  For  the 

4 x l014/cm  case  with  a 40  kA,  900  keV  beam,  we  have  that  the  initial 

growth  rate  is  given  by  6 = 3.8  X 10-3  x . Therefore,  from  Eq . (17)  we 

obtain  an  energy  delivered  in  the  initial  stages  of  the  instability  of 

As  = 5 X lO1"'  eV/cnT.  At  this  point  we  must  consider  the  effect  of 
o 

substantial  temperature  changes  due  to  heating  and  the  increasing  plasma 
density  due  to  ionization  of  the  3 X 1015/cm3  neutrals  present  in  the 
system.  If  we  combine  Eqs.  (15)  and  (17),  we  obtain  a wave  energy  depo- 
sition rate  appropriate  for  this  case 

ds  2 x 1023  kT  (t) 

— 1 = —st  • (18' 

dt  (n  (t)/4  x 1014  cm-3) 

P 

(in  this  equation,  kTg  is  an  effective  temperature,  since  there  may  be 
tails,  as  shown  in  Fig.  20.)  In  Fig.  15  we  saw  that  the  time  for  signi- 
ficant increase  in  density  in  this  case  is  comparable  to  the  pulse  dura- 
tion for  initial  energy  deposition  rates  of  order  1024  eV/cm3  - sec. 
Therefore,  it  is  reasonable  to  assume  small  density  change  and  take 
n ( t)  as  5 X 101'!/cm  . We  can  also  approximate  kTg(t)  by  10  eV.  This 
gives  an  "average"  value  for  the  wave  energy  deposition  rate  of 


de 

— - r*  1.5  X 1024  eV/cmJ 

dt 


sec  , 


a result  consistent  with  our  use  of  Fig.  15  for  guidance.  Thus,  the 
energy  deposition  via  waves  during  the  remaining  beam  pulse  is  Ae  « 

7.5  X 10l0  eV/cm2.  Using  Eq . (16)  to  obtain  a resistive  energy  deposition 
rate  equivalent  to  Eq . (18),  we  obtain 


de  1.6  x 101"  j 2 ( A / cm2 ) 

dt  ( n /b  X 1014  cm  3) 
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? X 
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eV 
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(19) 


the  latter  number  being  an  average  rate  during  the  pulse.  The  resis- 

tively  deposited  energy,  therefore,  totals  Ae  ~2.p  x 101'  eV/cnr. 

2 

Summing  the  three  components,  the  energy  deposition  according  to  the 
present  model  comes  to  1.0  X 101  eV/cnr.  Averaging  this  over  the  full 
beam  pulse  width  (70  ns),  we  obtain  an  average  energy  deposition  rate  of 
1.5  X 1024  eV/cm3  - sec,  greater  than  the  value  estimated  for  this 
case  from  Fig.  Ip  (including  the  40  eV  ionization  energy  and  the  2 V 2 
area  factor)  in  Sec.  IVB . 

Since  the  plasma  conditions  in  the  partially  ionized  cases  allow 
the  possibility  of  ion  acoustic  t -bulence,  we  must  consider  the  con- 
tribution which  this  instability  may  make  to  the  heating.  From  the  work 
of  Zavoisky,  et  al . , 7~  it  is  possible  to  obtain  a value  of  the  anomalous 

resistivity  due  to  ion  acoustic  turbulence,  ",  , for  v > v,  > c , where 

ia  e d s 

cg  is  the  ion  sound  speed,  (k(T  + T^)/M)  ^ . This  value,  in  a form 
appropriate  for  comparison  with  the  coefficient  of  js  in  Eq.  ( 1?) , is 
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„ _ l.S  X IQ17 eV 

'*‘a  (n^A  x 10* 4 cm  cm3  _ sec 
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Thus,  it  is  virtually  the  same  as  the  resistivity  due  to  the  e-e  mode, 
implying  negligible  additional  heating. 

The  final  case  to  be  discussed  here  is  the  one  with  an  initial 
density  of  8 x 1013/cm'  and  temperature  2 eV.  The  beam  was  1 MeV  and 
50  kA . At  this  low  plasma  density,  the  growth  rate  5 exceeds  the  colli- 
sion frequency  very  early  in  the  rise  time  of  the  current  pulse.  There- 
fore, we  take  half  the  beam  maximum  density  (5  x 1011/cm')  and  y = 2 in 

order  tc  obtain  the  6 to  use  in  calculating  As  . We  obtain  5/x  =s2  x 10 

o p 

and  Ae  ~ 1.0  x 10 1 eV/cirr  . During  the  next  stage  of  the  heating 

-3 

5/u.p  « 3 x 10  . However,  we  now  have  the  additional  complication  of  a 

significant  change  in  the  plasma  density  and  temperature  during  the 
course  of  the  interaction.  The  wave  and  resistive  energy  deposition 
between  20  and  70  ns  can  be  written 
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From  Fig.  A,  we  take  kT^  « 30  eV  as  an  average  over  that  time  interval, 
and  approximate  the  density  by  n^(t)  8 x 1 _1  + 0.12  t where  t 

is  in  ns.  This  gives 


55 


r 


Ae  X 1&'  eV/cm; ; Ae  « 4 x 101  ' eV  / cm3 , 

x 2 

for  a total  predicted  energy  deposition  of  3-5  X 1017  eV/crrr  . Assuming 
uniform  energy  deposition  over  the  20  cm^  beam  area,  this  is  a factor  of 
~ 3 greater  than  the  value  we  would  estimate  from  the  results  shown  in 
Fig.  14  (including  the  40  eV  ionization  energy  and  the  2 -/ 2 area  factor). 

Again  checking  the  possibility  of  ion  acoustic  resistive  heating, 
the  coefficient  of  jp  in  Eq.  (21b)  is  a factor  of  5 greater  than  the 
value  of  T.  at  n = c x 101'  /cm7 , obtained  from  Eq . (20).  However, 

71^a  “ 1/n  w^ereas  e-e  mode  resistivity  is  evidently  proportional  to 
n^  ‘ . Thus,  these  resistivities  are  about  equal  at  4 y 1014/cm  , which 

is  reached  at  about  50  ns.  Since  most  of  the  contribution  to  Ae  is  made 

2 

by  this  time (because  of  the  inverse  dependence  on  density  for  both  re- 
sistivities) , and  since  Ae  « Ae  , ion  acoustic  turbulence  will  again 
contribute  only  a small  amount  to  the  total  deposited  energy. 

In  Fig.  25,  we  plot  the  theoretically  predicted  energy  deposition 
as  a function  of  beam  to  plasma  density  ratio  for  the  three  cases  con- 
sidered. Also  shown  on  the  graph  are  the  experimental  results  from 
scattering  and  from  magnetic  diagnostics.  The  quantitative  agreement 
is  quite  reasonable  considering  the  difficulties  in  making  the  theo- 
retical estimates.  We  emphasize  that  this  comparison  is  a result  of 
including  the  early  and  late  stages  of  the  interaction,  and  is  not 
scaled.  By  contrast,  the  quasihydrodynamic  formulation  of  Thode  “ gives 
good  agreement  with  the  scaling  of  plasma  heating  by  an  IREB  as  a func- 
tion of  beam  and  plasma  parameters.  Quantitatively,  Thode  predicts  much 
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greater  energy  deposition  than  most  experiments  observe  (including  the 
present  one),  presumably  because  he  considers  only  the  early  stages  of 
the  interaction,  when  it  is  strongest. 

VI.  CONCLUSIONS 

In  this  paper,  we  have  presented  experimental  results  for  plasma 
heating  by  an  intense  relativistic  electron  beam.  The  initial  plasma 
density  ranged  from  ~ 5 X 10lj/ctr  to  ~ 3 x 10lo/cm"'  , the  lower  density 
cases  being  partially  ionized  and  the  higher  density  cases  highly 
ionized.  In  all  cases,  the  energy  coupled  from  the  beam  to  the  plasma 
is  greater  than  can  be  explained  by  binary  collisions  between  beam  elec- 
trons and  plasma  particles.  At  the  highest  plasma  densities,  classical 
damping  of  the  beam  induced  return  current  is  adequate  to  explain  the 
observed  heating.  However,  over  most  of  the  plasma  density  range  tested, 
i.e.  ^ 1.5  X 10lj/cm  , using  a 2 kA/cnr  beam,  the  plasma  heating  by  the 
beam  cannot  be  explained  by  classical  processes.  These  results  are  found 
to  be  explained  qualitatively  and  quantitatively  by  the  use  of  a full 
nonlinear  treatment  of  the  electron-electron  two  stream  instability  in 
the  kinetic  regime. 

The  results  of  this  experiment,  namely  2 - 7^/m  energy  coupling 

efficiencv,  at  beam  to  plasma  density  ratios  n,  / n in  the  range  10 

b p 

10  are  consistent  with  results  of  previous  experiments  with  comparable 
beam  and  plasma  parameters.3  ~ ’ Similar  conclusions  as  to  the 

energy  coupling  mechanism,  namely  the  electron-electron  instability,  have 
been  reached  in  most  of  these  experiments.  However,  the  use  of  more 
detailed  diagnostics,  and  a more  uniform  interaction  geometry  in  the 
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present  experiment,  and  a complete  nonlinear  treatment  of  the  e-e  mode 


instability  for  comparison  with  theory,  have  provided  strong  confirmation 
for  the  presence  of  this  instability. 

We  should  also  note  that  the  theory  we  have  used  predicts  that 
return  current  heating  due  to  the  parametrically  generated  ion  waves  W 
will  be  dominant  for  higher  beam  current  densities,  a fact  confirmed 
in  recently  presented  preliminary  results. 

Presuming  that  the  agreement  between  theory  and  experiment  is  not 
fortuitous,  it  is  interesting  to  note  a particular  implication  of  our 
theoretical  discussion.  Equation  (15)  for  the  wave  energy  transfer  rate 
during  the  stages  of  the  interaction,  when  combined  with  Eq.  (17), 
implies  a rate  which  increases  linearly  with  plasma  temperature  at  con- 
stant densities): 


de 

2_ 

dt 


52  — n kT 
P e 
ou 

P 


Thus,  as  the  plasma  heats,  the  coupling  efficiency  via  this  mechanism 
increases.  This  implies  that  the  use  of  1 us  pulse  IREB's  should  be 
advantageous  for  heating  plasma  in  a long  solenoidal  system  of  the  type 
suggested  by  Benford,  et  al.  Clearly,  this  equation  also  implies  that 
the  coupling  efficiency  of  the  entire  beam  will  be  increased  if  the 
plasma  can  be  preheated  to  a higher  temperature. 


?o 
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APPENDIX  A 


The  observation  of  the  scattered  power  spectrum  from  a beam  of 
laser  light  incident  on  a plasma  is  a well  established  technique  for  ob- 
taining information  on  the  electron  velocity  distribution  function. 

If  the  electrons  have  a nonrelativistic  (electron  temperature  « 1 ke' 
Maxwellian  velocity  distribution  function  with  a mean  velocity  of  zero, 
then  the  electron  temperature  Tg  may  be  obtained  by  plotting  the  loga- 
rithm of  the  measured  scattered  power  1^  vs  AX^  (X  - X ^ , where  X 
is  the  wavelength  of  the  incident  radiation  and  X ^ is  the  central  wave- 
length of  the  i1'1  channel.  Assuming  that  the  scattering  parameter 

a.  = > /2ttX„  « 1,  where  X„  = v fw  is  the  plasma  Debye  length  (v  is  the 
o D D e p e 

electron  thermal  velocity  and  id  the  plasma  frequency),  then  we  may  writ 


I . = I exp 

1 o 


c(a x)? 


where  I is  the  scattered  power  at  AX  = 0,  and 


4 v2V  sinp(0/2) 


■j.f  y 10' 

Xr  sin  (0/.°)  T (eV) 

o e 


The  velocity  of  light  is  c and  P is  the  angle  between  the  incident  wave- 

vector  and  the  scattered  wavevector. 

If  the  Maxwellian  distribution  has  a nonzero  mean  velocity  v , 

m 

then  the  scattered  power  spectrum  will  peak  at  a wavelength  Xc  such  thai 

(X  - X ) = (v  /c)X  • Then  instead  of  Eq . (Al1),  we  will  have 
' c o m o 


w 


I. 

1 


I ' exp 
o r 


(A3) 


where  now  Iq  ' is  the  scattered  power  at  X.  - X = 0,  and  G is  the  same 

so  long  as  I X - X I « X • 

1 o c o 

Suppose  now  we  have  five  scattering  channels,  labeled  1 thru  5, 
observing  scattered  power.  Then  the  ratios  of  the  observed  intensities 
will  be 


Taking  the  natural  logarithm  of  both  sides  of  Eqs.  (a4)  and  (A5),  solv- 
ing each  for  G and  equating,  we  can  solve  for  \ : 


X 


(X2  - X2)  An  - (x2  - X2)  An  ^ 

5 l 1 4 I 



I I 

2(X  - X ) An  — - 2(x  - X ) In  -S 
5 I 1 4 I 


(A5) 


Substituting  Eq . (Ao)  for  X back  into  either  Eq.  (a4)  or  (A5)  allows  a 

determination  of  G.  More  specifically  since  > and  0 are  known  and 

o 

t ^r\  I ° O 

fixed  ( o943  A and  90  , respectively,  in  the  present  experiment)  we  can 

solve  for  T : 
e 
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(A 


2.6  x 10'°  (X  + X - X - X ) (X  - X ) (X  - X ) 

■j.  1 4 g c;  i c i 

e II 

( X - X ) Xn  — - ( X - X ) Xn  — 

1 4 i 2 5 I 

5 4 

The  error  associated  with  the  use  of  X^  in  the  denominator  of  G in  place 

o r 

of  X , Is  small  relative  to  the  statistical  variations  of  channel  data 
c 

so  long  as  v « c. 

° m 

Once  X and  T are  known,  I , may  be  calculated  from  Eq.  (A3^  • 
c e o J 

Therefore,  the  electron  density  can  also  be  obtained  in  the  usual  way 
using  a Rayleigh  density  calibration. 

In  order  to  minimize  errors  in  Eqs.  (Ao)  and  (A7)»  we  clearly 
want  easily  measurable  intensities,  but  we  also  do  not  want  I « I 
and  I ss  I since  then  the  logarithm  of  the  required  ratios  will  be  neai 

O r~ 

c 

zero.  Placing  our  (equal  width)  channels  symmetrica lly  about  line 
center,  i.e.,  X -X  « - (X  - X ),  (X  -X)Ri-(X  - X ),  and  not 
using  I /l  and  I / I allows  satisfaction  of  both  of  these  requirements. 
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APPENDIX  B 


In  order  to  understand  the  observed  time  dependence  of  the  plasma 
density  and  temperature  in  the  partially  ionized  helium  experiments  de- 
scribed in  Sec.  IV7. B,  we  have  solved  the  following  set  of  coupled  differ- 
ential equations: 


dn 
o 

dt 


n n S + n n . ( Ct"*"  + of^n  ) , 
p o i p+  r 3 e 


(Bl) 


dn 


+ 


dt 


dn 
o 

dt 


|,-f-  | 1 

n n.S  + n n.,(a  + OL 

P + 2 p ++  r 3 


(B2) 


n-H-  = N ‘ 


(B5) 


n 

P 


= n++2 


(B4) 


_d_ 

dt 


- n kT 


= Q - 


n n S E 

P o i x 


n n.S  E - n n V] 
p + 2 2 p o ^ 


He- 1 


X.E. 

l l 


n n V X.E.  - P - Q - 
P + ^ J J B x+ 

He-II 


3 

2 


n kT  (n,a+ 
p e + r 


+ n 


h r 


) 


+ E 


+ n 


(B5) 


_d_ 

dt 


( 


2 n kT 

2 


+ 


) 


(Bo) 
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i (?■-"-) 


= Qj 


(B7) 


In  these  equations  we  are  considering  atomic  processes  only, 
neglecting  all  plasma  dynamical,  kinetic  and  transport  processes  (such 
as  expansion,  end  loss,  etc.).  This  is  equivalent  to  considering  the 
plasma  to  be  spatially  uniform  on  the  timescale  and  spatial  scale  of 
interest.  These  ideas  are  consistent  with  the  low  temperature  observed 
and  classical  transport  and  thermal  expansion.  In  solving  these  equa- 
tions we  assume  a Maxwellian  velocity  distribution  for  the  free  electrons 
and  ions. 

The  first  three  equations  describe  the  time  development  of  the 
neutral,  singly  and  doubly  ionized  helium  densities,  n , n and  n , 

O “T  "T+ 

respectively,  in  terms  of  the  electron  density  n^  and  the  various  rate 

coefficients.  S and  S are  the  ionization  rate  coefficients  for  neutral 
1 2 

+ ++ 

and  singly  ionized  helium,  respectively,  and  0:^  are  the  radiative 

recombination  coefficients  for  singly  and  twice  ionized  helium,  and  of*" 

3 

-f  | 

and  a are  the  three  body  recombination  coefficients.  Expressions  for 

3 

S and  S are  from  Drawin, /0  which  agree  well  with  the  calculations  of 
1 2 

Lotz.71  Adjustments  are  made  to  these  coefficients  to  account  for 
ionization  from  the  upper  levels  of  the  resonance  lines.  The  radiative 
recombination  coefficients  used  are  those  given  by  Seaton  ° (a  ) and 
Burgess  and  Seaton  1 (of^”),  and  a hydrogenic  expression  ' has  been 
utilized  for  the  three  body  recombination  coefficient.  Equations  (B3) 
and  (b4)  indicate  the  conservation  of  heavy  particles  (100  mTorr  He 
implies  N «3«5  X lO^/cm3)  and  charge  neutrality,  respectively. 
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The  last  three  of  the  above  equations,  when  combined  with  the 
other  four,  give  the  time  dependence  of  the  electron  temperature  T , 
and  the  singly  and  doubly  ionized  helium  temperatures,  T and  T , , , re- 

*i  "it" 

spectively.  The  terms  on  the  right  hand  side  of  Eq . (B5)  are,  in  order: 

1.)  the  heat  source  term  for  the  electrons,  usually  taken  to  be  a single 

resistive  term,  Tj2,  where  T is  the  resistivity  and  j is  the  plasma 

current  density;  power  loss  due  to  ionization  of  2.)  helium  atoms  and 

3.)  singly  ionized  helium,  where  the  ionization  energies  are  E =2 4.6  eV 

and  E = 54.4  eV ; power  loss  due  to  inelastic  electron  impact  excitation 
2 

of  4.)  helium  atoms  (He-l),  including  radiative  and  metastable  states, 

and  5.)  singly  ionized  helium  (He-Il),  where  E^  is  the  threshold  energy 

for  excitation  of  the  i*-^1  state  with  an  excitation  rate  coefficient  X^; 

6.)  power  loss  due  to  bremsstrahlung;  energy  transfer  rate  between  the 

electrons  and  J.)  singly  charged  ions  and  fc.)  doubly  charged  ions  (loss 

rates  so  long  as  T > T , T );  9*)  power  loss  due  to  radiative  recombina- 
e r "Hr 

tion  of  singly  and  doubly  ionized  helium;  and  10.)  power  gain  due  to  three 
body  recombination  of  singly  and  doubly  ionized  helium  where  E is  an 

average  energy  per  recombination  (taken  to  be  1/4  E and  1/4  E for  the 

1 2 

two  species).  The  excitation  rate  coefficient  X.  are  from  Drawin70  for 

l 

energy  loss  to  allowed  transitions.  The  coefficients  for  the  metastable 
states  of  helium  are  obtained83  from  measured  cross  sections.  The 

bremsstrahlung  power  coefficient  used  is  the  value  given  by  Spitzer 
and  the  electron  ion  energy  transfer  rates  are  obtained  from  the  classi- 
cal formula  given  by  Braginskii . 37 

The  resistivity  used  in  the  energy  source  term  Q = T|j2  in  Eq . (B5) 

consisted  of  the  sum  of  three  terms:  1.)  classical  resistivity  ’r 

J cZ 
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2.)  resistivity  due  to  ion  acoustic  turbulence,  71.  , as  determined  from 

la 

the  experimental  results  ot  Zavoiskiy,  etal./  and  3.)  resistivity  due 
to  electron-electron  two  stream  instability  turbulence,  77  , as  dis- 

cussed  in  Sec.  V.  Specific  values  used,  in  units  such  that  when  multi- 
plied by  j2,  with  j in  kA/cm2,  the  product  is  in  eV/cm  - sec,  were 


cl 


o x 


1022 


( 


j ?n  A 


.005  n 


(Be) 


ia 


5 y 1C 

n 


(B9) 


ee 


= 5 x lO- 


PS 


3/2 


(bio) 


— rj 

In  these  equations,  temperatures  are  in  eV  and  densities  in  cm  . The 
first  term  in  17^  is  from  Spitzer,f:'  In  A being  the  usual  coulomb  log- 
arithm, and  the  second  term  is  the  resistivity  due  to  electron- neutral 
collisions,  and  it  is  approximated  using  the  elastic  collision  cross 
section  given  by  Brown. a The  effective  collision  frequency  due  to  ion 
acoustic  turbulence  implied  by  Eq . (b9)  is  ~ ud^/BO,  where  uu  is  the 
plasma  frequency.  Since  Tg  » T+,  whenever  plasma  current  is  flowing 

in  the  partially  ionized  He  calculations  described  here,  T is  cut  off 

ia 

when  the  plasma  electron  drift  speed  drops  below  the  sound  speed 
c = JkT/H,  where  M is  the  mass  of  the  helium  ion,  as  per  the  dis- 

S 1 G 

cussions  in  Sec,  II. A and  V. 


6b 


By  far  the  most  Important  ionization  and  energy  loss  terms  in  Eqs . 
(Bl),  (B2)  and  (B5)  for  the  partially  ionized  helium  calculations  for 
which  results  were  given  in  Sec.IV.B  were  those  involving  ionization  of, 
and  line  radiation  from,  neutral  helium.  Therefore,  we  also  give  the 
specific  coefficients  used  for  these  terms,  in  units  of  eV-cm  /sec: 


S1  =2.35  X 10"7  Y (24.6)  exp  ^ - 24.6/Tg 


/f 


( B 11) 


v 

L 

He 


r X.E.  = 4.3  X 10"  ’ 1 0.276  Y (21.1)  exp  J^-  21. 


1/T 


(B12) 


0.2  Y (22.9)  exp 


- 22.0/Te  ] / nh  2.33  exp 


20/T 


(4.7  X 10"9  + 4.6  X IO"10  T - 1.2  x 10'11  T2)  , 

e e 


where  T is  in  eV  and 
e 

Y(X)  = ( + An  T 1 -2  (l  + T /X)l 

1 + X/T  \ 20  + X/T  *-  6 J 

e e 

The  first  set  of  terms  in  Eq.  (Bl2)  are  the  most  important  optical 
transitions,  and  the  second  set  of  terms  is  for  metastable  transitions. 

Equations  (Bl  - B7)  were  solved  using  a computational  scheme  de- 
veloped by  Boris  and  Winsor.90 

Calculated  electron  density  and  temperature  even  at  early  time 
(~  10  ns)  were  found  to  be  independent  of  initial  electron  temperature 
over  the  range  1 - 5 eV  and  ion  temperatures  from  0.1  - 1 eV.  Changes 
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of  the  order  of  25$  in  initial  density  did  affect  early  time  temperature 
(lower  initial  density  giving  higher  peak  temperatures  - see  Figs.  14 
and  15  - and  conversely),  but  only  small  changes  in  electron  density  and 
temperature  for  t > 100  ns  were  obtained. 
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Fig.  1 — A sample  set  of  oscillographs  of  (a)  the  Marx  generator  charging 
the  Blumlein  pulse  forming  line,  (b)  the  electron  beam  diode  voltage,  and 
(c)  the  electron  beam  current  in  the  diode 


icmatic  of  the  experiment  showing  the  location 
of  key  components  and  diagnostics 


(a)  PARTIALLY  IONIZED  CASE 
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(b)  HIGHLY  IONIZED  CASE 

Fig.  4 — Sample  oscillographs  from  a typical  channel  of  the  Thomson  scattering 
system  for  the  (a)  partially  ionized  and  (b)  highly  ionized  cases.  The  laser  timing 
r is  from  the  start  of  the  beam  induced  enhanced  continuum  to  the  peak  of  the 
laser  pulse.  The  large  enhanced  continuum  during  the  beam  pulse  in  the  par- 
tially ionized  case  required  r ^ 70  ns. 
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Fig.  6 — Determination  of  the  ion  temperature  from  the  Doppler 
broadening  of  the  HQ  line  (6531  A) 
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Fig.  7 — Plasma  density  versus  temperature  (from  Thor 
of  beam  shots  corresponding  to  the  plasma  only  shots  si 
dicate  two  data  points  with  unusually  high  temperatu 
emphasize  the  trend. 
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Fig.  8(a)  — Plasma  temperature  as  a function  of  time  (from  Thomson  scattering) 
for  shots  having  plasma  density  in  the  ranges  (a)  (7  ± 3)  X 1014/cm3,  (b)  (1.5  ± 
0.5)  X 1015/cm3,  and  (c)  (3.1  ± 1.0)  X 1015/cm3.  Parts  (a)  and  (b)  include 
data  points  (shown  with  error  bars)  obtained  during  a previously  published  run.25 
Arrows  on  error  bars  indicate  the  presence  of  tails  containing  an  unknown  total 
energy,  but  estimated  to  raise  the  average  electron  energy  to  near  the  top  of  the 
arrow.  The  curves  in  (c)  are  obtained  assuming  classical  resistive  heating  and  the 
plasma  current  densities  shown. 
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Fig.  8(c)  — Plasma  temperature  as  a function  of  time  (from  Thomson  scattering) 
for  shots  having  plasma  density  in  the  ranges  (a)  (7  ± 3)  X 1014/cm3,  (b)  (1.5  ± 
0.5)  X 1015/cm3,  and  (c)  (3.1  ± 1.0)  X 1015/cm3.  Parts  (a)  and  (b)  include 
data  points  (shown  with  error  bars)  obtained  during  a previously  published  run.25 
Arrows  on  error  bars  indicate  the  presence  of  tails  containing  an  unknown  total 
energy,  but  estimated  to  raise  the  average  electron  energy  to  near  the  top  of  the 
arrow.  The  curves  in  (c)  are  obtained  assuming  classical  resistive  heating  and  the 
plasma  current  densities  shown. 
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LASER  DUMP 


Fig.  9(a)  — Simplified  interaction  geometry  at  the  laser  scattering  part 
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Fig.  9(b)  — Intensity  versus  the  square  of  the  wavelength  shift  data  for 
selected  prep  shot  and  beam  shot  cases  to  illustrate  the  observed  shift 
from  blue  side  enhancement  to  red  side  enhancement 
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Fig.  10  — Calculated  scattering  line  center  (Appendix  A)  versus  (a)  density,  and 
(b)  temperature,  for  beam  and  plasma  shots.  Arrows  indicate  the  line  center  of 
the  17.1  eV  and  12.7  eV  cases.  The  tendency  for  beam  shots  to  have  a larger 
center  wavelength  is  evident  in  both  figures. 
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Fig.  10  (Continued)  — Calculated  scattering  line  center  (Appendix  A) 
versus  (a)  density,  and  (b)  temperature,  for  beam  and  plasma  shots. 
Arrows  indicate  the  line  center  of  the  17.1  eV  and  12.7  eV  cases.  The 
tendency  for  beam  shots  to  have  a larger  center  wavelength  is  evident 
in  both  figures. 
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Fig.  11(a)— Scattering  line  center  versus  (a)  density,  and  (b)  temperature, 
when  only  the  upper  half  of  the  scattering  volume  was  viewed.  Den- 
sity and  temperature  for  the  lower  half  of  the  scattering  volume  are 
shown  in  (c)  and  (d)  respectively.  The  upper  half  beam  results  show 
random  line  center  shifts  relative  to  prep  shot  results,  whereas  in  the 
lower  half  case  the  beam  shot  line  centers  are  shifted  to  longer  wave- 
length relative  to  prep  shot  centers. 
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Fig.  11(b)  (Continued)  — Scattering  line  center  versus  (a)  density,  and  (b)  tem- 
perature, when  only  the  upper  half  of  the  scattering  volume  was  viewed.  Den- 
sity and  temperature  for  the  lower  half  of  the  scattering  volume  are  shown  in  (c) 
and  (d)  respectively.  The  upper  half  beam  results  show  random  line  center  shifts 
relative  to  prep  shot  results,  whereas  in  the  lower  half  case  the  beam  shot  line 
centers  are  shifted  to  longer  wavelength  relative  to  prep  shot  centers. 
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Fig.  11(d)  (Continued)  — Scattering  line  center  versus  (a)  density,  and  (b)  tem- 
perature, when  only  the  upper  half  of  the  scattering  volume  was  viewed.  Den- 
sity and  temperature  for  the  lower  half  of  the  scattering  volume  are  shown  in  (c) 
and  (d)  respectively.  The  upper  half  beam  results  show  random  line  center  shifts 
relative  to  prep  shot  results,  whereas  in  the  lower  half  case  the  beam  shot  line 
centers  are  shifted  to  longer  wavelength  relative  to  prep  shot  centers. 


92 


r 


time  (nsec) 

Fig.  12(b)  (Continued)  — Typical  Faraday  cup  oscillogram  for 
the  magnetic  probe  experiment 


Fig,  13(a)  — The  change  in  the  axial  magnetic  field,  SBZ  as  a function 
of  probe  position  at  various  times  during  the  beam-plasma  interaction. 
Shot  to  shot  reproducibility  is  illustrated  by  individual  data  points 
shown  at  one  radial  position  in  each  figure. 
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Fig.  13(b)  (Continued)  — Horizontal  magnetic  field,  Bx,  as  a function 
of  probe  position  at  various  times  during  the  beam-plasma  interaction. 
Shot  to  shot  reproducibility  is  illustrated  by  individual  data  points 
shown  at  one  radial  position  in  each  figure. 
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Fig.  14(a)  — Results  for  the  8 X 1013  case  showing  plasma  electron  temperature, 
data  as  a function  of  time.  Typical  error  bars  are  shown.  Solid  curves  are  ob- 
tained from  the  numerical  model. 


Fig.  14(b)  (Continued)  — Results  for  the  8 X 1013  case  showing  density, 
data  as  a function  of  time.  Typical  error  bars  are  shown.  Solid  curves  are 
obtained  from  the  numerical  model. 
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Fig.  15(a)  — Results  for  the  4 X 1014  case  showing  plasma  electron  temperature, 
data  as  a function  of  time.  Typical  error  bars  are  as  shown.  Solid  curves  are  ob- 
tained from  the  numerical  model. 
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Fig.  15(b)  (Continued)  — Results  for  the  4 X 1014  case  showing  density, 
data  as  a function  of  time.  Typical  error  bars  are  as  shown.  Solid  curves 
are  obtained  from  the  numerical  model. 
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Fig.  16  — Typical  corrected  diamagnetic  loop  signal  for  the  8 X 1013  case 


101 


COMPENSATING 

COIL 


TIME  (nsec)  — - 

(b) 

Fig.  17  — (a)  Circuit  used  to  approximate  the  diamagnetic  loop  circuit, 
(b)  Expected  signal  from  the  diamagnetic  loop  resulting  from  the  trans- 
verse energy  pulse  described  in  the  text  and  using  the  circuit  shown  in 
part  (a). 
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Fig.  18  — Diamagnetic  loop  signals  from  loops  at  both  ends  of  the  theta  pinch 
in  a partially  ionized  hydrogen  experiment.  Oscillations  were  present  in  most 
hydrogen  cases. 


Fig.  19(a)  — Comparison  of  diamagnetic  loop  amplitude  and  beam  energy  on  the 
calorimeter,  for  partially  ionized  hydrogen  experiments 
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Fig.  21  — Beam  current  density  pulse  shape  used  for  the 
theoretical  comparison  with  experiment 


M 


10 


17 


ro 

6 

§ I016 

<u 


>- 

O 

cr 

UJ 


u 


Q|0I5 


LlI 


CO 

o 

0. 

UJ 

Q 


10 


14 


TIME  (nsec) 


Fig.  22  — Beam  energy  deposition  rate  as  a function  of  time  for  the  highly 
ionized  case.  Continuous  lines  are  averaged  results  from  the  approximate 
analytic  theory.  Line  (a)  represents  deposition  due  to  wave  damping,  while 
(b)  is  the  total  including  the  resistive  heating.  The  points  are  results  from 
the  numerical  solution  of  the  exact  mode  coupling  equations  including  only 
wave  damping. 
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Fig.  23  — Time  dependence  of  anomalous  resistivity  for  the  case  considered 
in  Fig.  22.  Again  the  line  is  from  the  analytic  model  while  the  points  are 
from  the  numerical  experiments  averaged  over  a ~ 1 ns  time  interval. 
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Fig.  25  — Energy  deposition  as  a function  of  plasma  density:  a comparison 
of  experimental  observations  and  theoretical  predictions.  Note  that  the  two 
lower  density  cases  are  partially  ionized  helium  and  the  other  three  are 
highly  ionized  hydrogen. 


